General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



o* w//j/yits 

it 


(NASA-CR- 173924) THEORETIC AL .RESEARCH 
PROGRAM TO PREDICT THE PROPERTIES OF 
MOLECULES AND CLUSTERS CONTAINING TRANSITION 
METAL ATOMS Final Report, Jul. 1981 - J'ul. 
1984 (Polyatowics Research, Inc.) 97 p 


N84-32502 


Unclas 

01066 


Theoretical Research Program to Predict the 
Properties of Molecules and Clusters containing 
Transition Metal Atoms 


Final Report Covering Period July 1981 - July 1984 


Supported under NASA Cooperative Agreement NCC2-148 


NASA Technical Officer: Dr. David Cooper 

NASA Ames Research Center 
STC 230-3 

Moffett Field, CA 94035 

I 

Principal Investigator: Dr. Stephen Walch * 


Polyatomics Research Institute 
1101 San Antonio Road, Suite 420 
Mountain View, CA 94043 


The enclosed document is the final technical report on 
work performed under NASA Cooperative, Agreement NCC2-148. 
This report covers the work of Dr. Stephen Walch who acted 
as principal investigator. The information contained 
herein is complete and was submitted in July, 1984. 


TABLE OF CONTENTS 


INTRODUCTION 1 

TRANSITION METAL HYDRIDES . 2 

TRANSITION METAL DIMERS 6 

Fe H CLUSTER CALCULATIONS 15 

n 

LIST OF PUBLICATIONS 16 

APPENDIX A. . . * 19 


APPENDIX B - Reprints of papers describing calculations 
for Transition Metal Hydrides and Dimers. 


X. INTRODUCTION 


The primary focus of this research has been the theoretical study 
of transition metal (TM) chemistry, A major goal of this work is to 
provide reliable information about the interaction of H atoms with iron 
metal. This information is needed to understand the effect of H atoms 
on the processes of embrittlement and crack propagation in iron. The 
method employed in the iron hydrogen studies is the cluster method in 
which the bulk metal is modelled by a finite number of iron atoms. 

There are several difficulties in the application of this approach to 
the hydrogen iron system. First the nature of TM-TM and TM-H bonding 
for even diatomic molecules was not well understood when these studies 
were started. Secondly relatively large iron clusters are needed to 
provide reasonable results. Therefore it is not possible to include all 
of the valence electrons on each iron atom. In the present calculations 
only the 4s electron of each iron atom was included directly in the 
calculation while the Ar core and 3d electrons are incorporated into the 
effective core potential (ECP) based on the 4s 3d state of the Fe atom. 

In order to understand the effects of removing the 3d electrons from the 
valence space, it is necessary to have detailed information about the 
nature of TM-TM and TM-H bonding. For these reasons we have carried out 
in parallel two levels of calculation: i) highly accurate studies of TM 

dimers and hydrides. (These studies are carried out with large basis sets 
and include extensive electron correlation) . ii) studies of large clusters 
of Fe atoms with and without H adatoms. [These studies are carried out at 
the SCF level and make the approximation of Incorporating the TM core and 
3d electrons into the core (ECP approximation) . ] The studies of diatomic 
species provide benchmark calculations which are used to test the 
approximations in the cluster studies, while the cluster studies on Fe^H 
are designed to provide potential function input for studies of hydrogen 
embrittlement. 

Since much of the work to be discussed here has been described in 
manuscripts, relevant manuscripts have been included in the appendices 
and the reader is referred to these papers and to the publications list 
for further details. This report concentrates on the TM hydrides, TM 


dimers, and Fe H cluster results. The Fe H cluster results are 
n n 

described in detail in a manuscript included as Appendix A. The 
results for the TM hydrides and dimers are summarized in Sections IX 
and III, respectively; while, Appendix B contains copies of publications 
on these topics. During the course of this work several other projects were 
undertaken. These include all-electron calculations on the CsH and CS2 
molecules, calculations on the alkali dimers, an ab-initio study of core- 
valence correlation, and studies of atomic correlation and basis sets. 
References to this work are given in the publications list. 


II. TRANSITION METAL HYDRIDES 


2 n 1 n*fl 

Fig. 1 shows the relative ordering of the 4s 3d and 4s 3d 

states of the TM atoms. The Important features of Fig. 1 are: 1) a 

1 n+1 

monatonlc decrease In excitation energy from Sc to Cr with 4s 3d n an 

excitated state for Sc through V but the ground state for Cr. 11) a 

sharp reversal in the ordering of the states for Mn with 4s 3d n again 

lower, and ill) a monatonlc decrease in excitation energy from Mn to 
1 10 

Cu with 4s 3d the ground state of Cu. The character of the ground 
states of the TM hydrides is strongly correlated with the ordering of 

atomic states. For example, for CrH the ground state is ^E* arising 

15 7 + 2 5 

from 4s 3d while for MnH the ground state is E arising from 4s 3d . 

2 n ^ H"f*T 

For elements where the 4s 3d and 4s 3d states are closer in energy 

5 

a strong admixture of states occurs, e.g. for VH the ground state is A 

2 2 1 112 1 

arising from a mixture of 4s 3dTT 3d6 and 4s 3da 3d7T 3dt$ atomic character 
with a 3d population of A;3.4. 

The ratio of sizes of 4s and 3d orbitals Increases monatonlcally 

from 2.364 to 3.239 going from Sc to Cu. Because the 4s orbital is 

so much larger than the 3d orbital we find the bonding in the TM hydrides 

involves primarily the 4s electrons (with some admixture of 4p character) 

while the 3d electrons remain essentially atomic like. The one exception 

1 + 

is ScH which shows a Sc(3d)-H(ls) bond in the lowest E state. Note 
that Sc is most likely to show 3d bonding since the 4s and 3d orbital 
sizes are more nearly comparable. 

CASSCF/CI studies have been carried out for T1H, VH, CrH, MnH, 

FeH, and NiH. These calculations show good agreement with experiment for 
the R^, o) e » and D g ; thus confirming the theoretical model. 

Studies have also been carried out for the state of FeH - . 

These calculations show that photodetachment from FeH - leads to a large 

4 6 

change in R g for the A state but little change in R^ for the A state. 

These results confirm the assignment of Stevens, Fiegerle, and Lineberger 

(SFL) of the vertical photodetachment transition with a long vibrational 

4 4 

progression to A, and hence that the ground state of FeH. is A. However 

the calculations show that the simple theoretical model used by SFL to 

4 

assign the above transitions is not correct in that the A state of FeH 


and the 5 A state of FeH~ have 3d populations of 6.26 and 6.50 respectively, 
whereas SFL assumed 7.0 for both of these states. Thus the theoretical 
calculations In conjunction with the experiment by SFL establish that 

4 

the ground state of FeH is A. 
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-!“• 1* 4s^3d n+ ^-4s^3d n excitation energies of * 

scandium to copper [E(sd n+1 )-E(s 2 d n ) ] . All units 
are in eV. 



III. TRANSITION METAL DIMERS. 


For the TM dimers molecular states may arise from three different 
atomic asymptotes J 4s^3d n + 4s 2 3d n , 4s^3d n + 4s’*‘3d n ^"‘ 1 ', and 4s‘ L 3d n ^ 1 + 
4s^3d n ^. The accessibility of these asymptotes of course depends on 
the 4s^3d n+ ‘*‘ - 4s^3d n excitation energies (see Flg.l). Thus for SCg 
we find that the low-lying states arise from the first two asymptotes 
with 4s^3d n+ ^ + 4s^3d n+ ^ too high in energy to lead to the ground state 
while for Tig* Vg, and Crg the ground states arise from the 4s^3d n ^4-4s^3d n+ ^ 
atomic asymptote. 

Because the 4s orbital is significantly larger than the 3d 

orbital the predominant interaction in the TM dimers at large internuclear 

distance (R) is between the 4s electrons with very little 3d Interaction. 

2 n 2 n 

For states arising from the 4s 3d + 4s 3d atomic asymptote this 

interaction is basically repulsive and only a shallow well at large R g 

('v8.0a Q ) arising from the 4s^4p near degeneracy effect is observed. For 

states arising from the 4s^3d n + 4s'*'3d n ^ atomic asymptote the 4s 

interaction is weakly bonding at intermediate R but is repulsive at 

small R leading to intermediate R £ values (^5.0a Q ). At these R g values 

the 3d-3d overlaps (S) are small which favors one-electron 3d bonds 

whose bonding terms vary with distance like S over two-electron 3d bonds 

2 

whose bonding terms vary with distance like S . Finally for states arising 

1 n-fl 1 n+l 

from the 4s 3d + 4s 3d atomic asymptote the 4s-4s interaction is 
attractive and also appears to be relatively flat well Inside the optional 
4s-4s bonding radius. Thus, states arising from this atomic asymptote 
are able to move into short R^ regions where the 3d-3d overlaps are large 
enough to favor two-electron 3d bonding. 

For the TM elements with more than half filled 3d shells, the 
formation of 3d bonds becomes much less favorable for two reasons. First 
the 4s to 3d orbital sizes are larger for the right half of the first 
transition row. Secondly the presence of doubly occupied 3d orbitals 
leads to repulsive interactions which effectively cancel any bonding 
interactions from the 3d shell. Thus the bonding here is dominated by 
the 4s electrons. An example is Cug where the Cu atom has a 4s^3d^ 


ground state and the bonding Is predominantly a 4s-4s bond with the 
3d electrons remaining atomic like. 

The dominant configurations for important states of those TM 
dimers which have been studied are: 
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Fig. 2 shows calculated potential curves for low-lying states of 

Sc_. The initial study of Sc 9 by Walch and Bauschlicher found only weakly 

2 1 2 2 
bound states arising out of the D + D asymptote in contrast to mass 

spectrometric experiments which indicated strong bonding (D = 26± 5kcal/mole) . 

5 e 2 4 1 1 1 

A z state was found which was bound by 'v-O.SeV with respect to D + F(4s 4p 3d ) 

u 2 2 
atomic limit, but unbound with respect to D + D. However at about the same 

time that this work was published, matrix isolation studies by Knight, VanZee 

and Weltner indicated a bound ^ state of Sc„. From the ESR studies it 

2 1 4 .,1 2 

appeared that this state arose from the D + F (4s 3d ) atomic asymptote 
which had not been studied in detail in the previous theoretical studies. 

A re-examination of this system revealed a new state which had been missed 
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lri the previous study because Its R (<v5,0a ) Is much shorter than the 

e 5 - 


R values for the states studied previously ('v? *0» o ) . The E u state 


turned out to be of considerable theoretical interest because it 
exhibited multiple 3d bonding (three one-electron bonds) and constituted 
the first theoretical evidence of multiple 3d bonding in a first row 
transition metal dimer. 


The studies on SCg were extended to Tig, Vg and Cr g. Initially 


74 . 4 4 5- 7 + 

the high spin ' E^ state of Tig was studied. The SCg E u and Tig E u 


'2 n 1 n-fl 

states come from the mixed (4s 3d +4s 3d ) asymptote and exhibit 


multiple one-electron 3d bonding. These states have long bond lengths 

-1 5 - 

(R >5.0a ) and small vibrational frequencies (w o-200cm ). The E 

' e nj o c e a 

state of SCg is consistent with the E state which has been observed in matrix 
studies by Knight, Van Zee, and Weltner. The bond length here is not known 


-1 


but the calculated vibrational frequency of 184cm is reasonably close to 
the experimental value of 238. 9cu ■*-. 


For Ti 2 from reference to Fig 1 one sees that the excitation energy 
to the 4s 1 3d n *’ L +4s 1 3d n+1 atomic asymptote is about half as large as for 


- , 1 rt 1 O 

Sc„ . Thus the 1 E + state which arises from the 4s 3d J +4s 3d atomic asymptote 
• g 


becomes a competitor for the gound state of Tig. The bonding here is a 


triple two-electron 3d bond (3da» 3dir , 3dir ) . This leads to a short R 

-i x y 

state, R = 3.73a and w “ 438cm A . The bond length is not known 
e o e 


-1 


experimentally but the experimental vibrational frequency is 407.9cm 
which is consistent with the 1 E g state of Tig, but inconsistent with the 

^E + state. Based on this we tentatively assign the ground state of Tig 

U 1 + 

35 Jj * 

g 3 - 

Fig. 3 shows calculated potential curves for the E g state of Vg. 


O 

The e state of V„ has the same triple two-electron 3d bond as in Tig 

with the remaining two electrons in the 3dfi orbitals. Because the 3ds 

orbitals still have small overlaps in the region near R g (^3.5a Q ), the lowest 

state is a triplet state arising by forming two one-electron 3dfi bonds. The 

R and m values obtained from the CASSCF curves are in good agreement with 
e e 

the recent results of Langrldge-Smlth, Morse, Hansen, Smalley and Merer for 

O 

the J E ground state of Vg. An important feature of Fig. 3 is the large 
effect of 4f functions, an effect which is also evident in the Grg curves 
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FIG. 2. Fotential energy curves for the low-lying 
states of Sc2 from CASSC? CICSP) calculations. The 
locations of the 2 D + 2 D, 2 D + 4 F(As^4p^3d^) , and 
2 D + 4 F(4s^3d 2 ) asymptotes are indicated. 
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In Fig. 4. Since these states exhibit strong 3d bonding, the large 

effect (H.Oev) for 4f as a polarization function is not surprising. This 

5 - 

large effect is not observed for the one-electron 3d bonds in the Sc 9 £ 

7 i £ U 

or Tig states. However, this is expected given the strong R 

dependence of this effect. For example from Fig, 4 one sees that for 

Cr_ near R the effect of 4f is large but is near zero by 3.75a , 

Z e o 

Fig. 4 shows the calculated potential curves for Cr^. Here the 
3dd orbitals are doubly occupied which is expected to be unfavorable 
based on the V 2 result that the 3d« orbitals were preferentially singly 
occupied. This is consistent with the weaker bonding in Crg, & e M..0eV 
as compared to V^, D e ^2,5eV. Because of the weaker bonding in Cr ^ the 
CASSGF potential curve is not hound. However the potential curve does 
exhibit a shoulder near the experimental R £ which is suggestive of an 
inner well. 

A major problem with the CASSCF studies for Vg and Cr^ is that 
only a small percentage of the binding energy is obtained for V ^ and no 
well is obtained for Cr^. This result is not unexpected for CASSCF. 

Normally these problems would be corrected by configuration interaction 
(01) . However, the Cl expansions required for and Cr^ exceed current 
computational capabilities. Several different attempts have been made to 
include the missing correlation in other ways. One approach due to Goodgame 
and Goddard assumes that the missing correlation serves mainly to correct 
the location of the ionic atomic asymptotes. These authors attempt to include 
these effects by empirical modification of the integrals to correct the 
atomic ionization potentials and electron affinities to agree with experiment. 
This method does lead to a reasonable potential curve for Cr^ although the 
bond length is somewhat too short which suggests that this method over 
corrects. Another approach by Walch attempts to include extra correlation 
effects by expanding the valence space in the CASSCF calculation. Here 
atomic 4p and 3d' terms are added (where 3d' is a tight diffuse correlating 
orbital for the 3d). This approach should include the principle correlation 
effects needed to describe charge transfer within the 4s and 3d shell. 
Unfortunately, it is not possible to add all these extra valence orbitals at 
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once (due to computational limitations) and the extra orbitals are added 
separately by symmetry blocks. From checks on Hg it appears that there 
is no real problem with additivity due to this approach. However, cheeks 
on Ng and Tig indicate that because the method dissociates to SCF atoms 
certain atomic correlation terms are not included and these terms are 
more Important in the large R region than they are near R fi . Thus, the 
method does overestimate the binding energy that would be obtained in a 
more complete MCSCF calculation. In spite of these difficulties the 
method does nearly reproduce the Cl potential curves for Hg , Ng and Tig 
and therefore we reproduce the estimated potential curve for Crg in Fig. 5, 
Analysis of this potential curve gives (experimental values in parenthesis) 


R g = 1.78A°(1.68A°) , 


to = 383cm ^(480cm ^) and D = 0.71eV). Note that 
e e 


the long bond length and small to are consistent with underestimating the 

6 

binding energy. 

A significant feature of the studies of the TM dimers is the presence 

of outer wells associated with 4s-4s bonding and one-electron 3d bonding 

7 + 

in some cases (e.g. the I state of Ti„ and analogous states in V 9 ) , and 

u ^ 1 + Z 
inner wells associated with two-electron 3d bonding (e.g. the 2 state 


of Tig and the 


3„- 


£ state of V_) . These types of effects have been observed 
g 2 ' 

experimentally for small Fe clusters where the dimer exhibits a bond length 


of ^3.83^^ which presumably involves some 3d bonding but larger clusters 

exhibit longer bond lengths approaching the nearest neighbor distance in 

BCC Fe (^4. 7a ) . 

o 
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7IG. 5. Estimated potential curves for based on 
extended valence space CASSC7 calculations. The effect 
of 4p is estimated based on calculations on T^. 
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IV. Fe H CLUSTER CALCULATIONS 
n 

The Fe^H cluster calculations are described in detail in a 
manuscript contained in Appendix A. For a summary of the main results, 
the reader is referred to the introduction of the above paper. 
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APPENDIX A 


The following paper discusses the Fe n H cluster calculations. 


Model Studies of the Interaction of 
H Atoms with BCC Iron 


Stephen P. Walch 

4 * 

Polyatomics Research Institute 
Mountain View, CA 94043 

Abstract 

Ab Initio/Effective Core Potential cluster studies are reported for the 

interaction of H atoms with BCC iron. The calculations use a one-electron 

1 7 

ECP based on the 4s 3d state of the Fe atom. Two-fold and four-fold sites 

on the (100) surface as well as octahedral, tetrahedral, and trigonal interior 

sites were studied. Four-fold surface sites are found to be bound by 'vL.SeV 

with the H atom 'uO.Sa above the surface. Penetration of the surface at a 

o 

four-fold site involves movement toward a second layer atom and is expected to 
be unfavorable. Two-fold surface sites have small binding energies 'v0.25eV. 
Penetration of the surface at this site involves movement toward a tetrahedral 
interior site and is downhill in energy. Tetrahedral interior sites are found 
to’ be bound by o.l.3eV and are a minimum on the potential energy surface. 
Octahedral sites are a maximum on the potential energy surface and are estimated 
to be r v0.2eV higher (including lattice relaxation effects). Trigonal sites 
are found to be a saddle point connecting adjacent tetrahedral sites and this 
pathway leads to an estimated barrier to diffusion of 'vO.leV (including lattice 
relaxation effects) . The volume expansion for a H atom in a tetrahedral site 
is calculated to be 21%. 

* Supported by NASA grant // NCC2-148 
t Mailing address: 1101 San Antonio Road, Suite 420 


dimers* and Fe H cluster results. The Fe H cluster results are 
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described in detail in a manuscript included as Appendix A. The 
results for the TM hydrides and dimers are summarized in Sections II 
and III, respectively; while, Appendix B contains copies of publications 
on these topics. During the course of this work several other projects were 
undertaken. These Include all-electron calculations on the CsH and 
molecules, calculations on the alkali dimers, an ab-initio study of core- 
valence correlation, and studies of atomic correlation and basis sets. 
References to this work are given in the publications list. 
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X, Introduction 


Hydrogen embrittlement of metals is an important technological problem 
11]. As part of a program to understand the effect of H atoms on crack propa- 
gation in iron we have undertaken a theoretical study of the interaction of 
H atoms with clusters of iron atoms in the BCC structure (a-Iron) . In the 
studies reported here we first fixed the Fe atoms at the lattice positions 
of the perfect crystal and studied the interaction with a H atom at two-fold 
and four-fold sites on the (001) surface and at a series of octahedral, tetra- 
hedral, and trigonal sites interior to the solid. These studies determine 
the geometry and binding energy at surface chemisorption sites and provide 
information on the barriers to diffusion within the metal and to penetration 
of the surface at two-fold and four-fold sites. We then relaxed the geometry 
of the nearest neighbor Fe atoms for the octahedral, tetrahedral and trigonal 
interior sites. These studies are important because they further refine the 
energetics for the diffusion process and because significant expansion of the 
lattice is found for tetrahedral sites. This expansion could induce stress 
in the metal leading to fracture. 

The computational method employed here is the ab initio SCF method using 

a one-electron (As valence) iron effective core potential (ECP) [2] based on 

17 7 

the As 3d state of the Fe atom. Here the Ar core and the 3d configuration 

are incorporated into the ECP leaving only a single As valence electron. 

The justification for including only the As electrons in the valence space 

comes from all-electron studies of diatomic molecules Involving transition 

metal (TM) atoms. Here is is found that for elements on the right hand side 

of the first transition row, the bonding in the hydrides [3] and dimers at 

the nearest neighbor distances in the metal [A-6] Involves the As,Ap electrons 

with the 3d electrons remaining essentially atomic like. The selection of 

the As 3d atomic configuration of the Fe atom is based on experimental bulk 

magnetic information [7] and on calculations [8] which show that although the 

2 6 17 

As^3d° state is the ground state of the free Fe atom with the As A 3d state 

1 7 

0.88eV higher, in an environment appropriate to Fe metal the As 3d state 
becomes the lowest state. 


The present calculations Involve Fe clusters of up to 66 Fe atoms and 
show slow convergence with respect to cluster sizej cluster edge effects still 
introduce significant uncertainties even for the largest clusters studied. 
Important qualitative features of the H~Fe n interaction include: i) H atoms 

are most stable at surface four-fold sites (binding energy a/1.5eV) with a 
barrier to movement into the bulk (toward an Fe atom in the second layer), 
ii) Two-fold sites are only weakly bound but it is downhill from this location 
to move into the bulk (toward an Interior tetrahedral site), lii) Interior 
octahedral and tetrahedral sites are less stable than four-fold surface sites 
(binding energy n/l.SeV) with tetrahedral sites lower in energy (n*0.23eV repara- 
tion for the unrelaxed lattice and 'uO.lSeV separation for the relaxed lattice 
based on Fe^H) * The tetrahedral site is a minimum on the potential energy 
surface while the octahedral site is a maximum, iv) Tetrahedral Interior sites 
show a large volume expansion ^21%. v) The lowest pathway for T«-*-T diffusion 
involves a trigonal site which is a saddle point on the potential energy surface 
(Calculated barrier 'vO.neV for the unrelaxed lattice and 'VO.lOeV for the relaxed 
lattice based on Fe^H) . 

Section II discusses some features of the electronic structure of Fe 
clusters. Section III discusses the cluster calculations. Section IV dis- 
cusses the features of the calculated potential, while Section V compares 
the computed energetics to experiment. 


II. Electronic Structure of the Fe Clusters 


For the TM atoms of the first transition row the 4s orbital is signifi- 
cantly larger than the 3d orbital. As one moves from left to right in the 
first transition row the ratio <r 4 8 5 ^ <r 3d > increases monatonically from 
2.36 to 3.24 with a value of 2.95 for Fe [9], Thus, for large bond distances 
one expects bonding to the 4s orbital to be more favorable than bonding to 
the 3d orbital W3>th the 4s becoming increasingly favorable on the right side 
of the row. 

For the TM hydrides only ScH [10) shows a short R state involving 3d 

bonding; whereas the other hydrides including FeH show predominately 4s 

bonding character with the 3d orbitals remaining quite atomic like [3]. The 

situation for the dimers is similar with the possibllty for elements on the 

left half of the row of 4s-4s bonding at large R and 3d-3d bonding at small R. 

Examples of the latter are [11] and [11,12). With more than half 

filled 3d shells the favorability of 3d bonding decreases since the repulsion 

between doubly occupied 3d orbitals cancels 3d bonding interactions arising 

out of the singly occupied 3d orbitals. For Cu 9 [4] with a closed 3d shell 
1 10 * 

(for the 4s 3d state) the bonding is described as predominately 4s-4s bonding 
with the 3d orbitals essentially atomic like. Given this, one expects the 
dominant bonding interactions for Fe-Fe and Fe-H bonds to involve the 4s, 4p 
electrons of the Fe atom with the Fe 3d electrons remaining essentially atomic 
like. 

Fig. 1 shows the results of all-electron calculations for the Fe„ molecule. 

2 6 1 I 

The electronic state considered here arises from the 4s 3d + 4s 3d atomic 

limit. The significant feature of Fig. 1 is that there is weak 3d bonding 

in the small R region as evidenced by the low-spin Btate (maximum 3d bonding) 

being lower, but in the large R region the high spin state (no 3d bonding) 

is lower thus indicating negligible 3d bonding. Similar conclusions have 

l.J 


been reached by Shim and Gingerich [6] for mates derived from the 4s 3d *f 
4s *3d^ atomic limit. Thus we conclude that at the Fe-Fe distances involved 
in BCC Fe metal-metal 3d bonding is of negligible Importance. (The nearest 


neighbor distance is 4.68a Q in BCC Fe [13)). For these reasons we believe 
it is a reasonable though extreme approximation t& include only 4s electrons 
In the valence shell with the 3d electrons incorporated into the ECP. 


3 
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The ECP Is based on the 4s 3d Fe atomic configuration. The choice of 

this configuration is supported by experimental magnetic date [7] which 

show 2.22 effective Bohr magnetons per Fe atom. Assuming that this magnetic 

moment arises only from the 3d electrons and that the 3d electrons are completely 

1 7 

high spin coupled, this magnetic data is co sistent with a mixture of 4s 3d 

8 17 

and 3d . However this data could also be consistent with 4s 3d and some 

1 7 

low-spin coupling. The choice of 4s 3d is also supported by calculations 
in which one Fe atom is placed in the center of a cubic arrangement of eight 

one-electron 4s*3d' r like Fe atoms (BCC unit cell). The result is that the 

17 2 6 

4s 3d configuration is 0.25eV below the 4s 3d configuration in this environ- 

ment, although for the free atom 4s 1 3d / is 0.88eV above 4s z 3d, Given this 

choice of atomic configuration with the 3d electrons incorporated into the ECP 

leaves only one valence electron per iron atom. This approach is very similar 

1 9 

to the one-electron Ni cluster studies based on the 4s 3d Ni atomic configu- 
ration [14]. 

Comparison between ECP and all-electron results can be made for FeH. 

Here the ECP calculation leads to R = 2.96 a and D = 1.85eV whereas the 

e o e , 

all-electron calculation leads to R = 3.02 a and D = 3.20eV for the ^A 

e o e 

state of FeH [5). However, a direct comparison here is not appropriate since 
the ECP calculation corresponds to pure 4s*3d^ atomic character whereas the FeH 
molecule in the ^A state is a strong mixture of 4s*3d^ and 4s^3d^. The latter 

configuration mixing, which occurs in the all-electron calculations, is ex- 

2 6 

pected to lengthen R^ (due to more 4s 3d character) and increase D g as 

compared to the ECP calculation; these expectations are borne out by the results 

1 7 

quoted above. Based on these results, to the extent that mixing of 4s 3d 
2 6 

and 4s 3d is important for the cluster, we expect the calculations with the 
ECP to slightly underestimate FeH bond lengths and to significantly under- 
estimate H atom binding energies. Relative energetics for a given cluster 
would still be expected to be accurate and the emphasis here is therefore 
on comparisons of relative energies for different H atom locations within a 
given cluster. 
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Ill, The Cluster Calculations 

The basis set for Fe is a (4s3p) / [ 2slp ] basis. The 4s functions are a 

(31) contraction of the outer four Wachters' 4s functions [15] . The contrac- 

2 

tion here is based on an atomic SCF calculation for the S state arising from 
the 4s* configuration. Note that this contraction is different from the all 
electron contraction because the ECP 4s orbital does not have a nodal structure. 
The 4p function was obtained by optimization of a 3 term GTO fit to a 2p STO 
[16] for the P atomic state arising out of the 4s configuration. The re- 
sulting exponents were multiplied by 1.5 to make them more suitable for de- 
scribing 4s + 4p correlation [17]. The H basis set is a (5slp)/[3slp] basis. 

The s basis is a (311) contraction of Huzinaga's 5s set [18], while the 2p 
functions are a single set of gaussian primitives with exponent 1.0. The one- 
electron ECP is given in Table I. 

As has been well documented [19] 4p functions are very important for 
describing the bonding in Ni clusters and these effects are also found to be 
important for Fe clusters. Omission of the 4p functions leads to a loss of 
most of the metal metal bonding energy and leads to overestimation of metal-H 
binding energies due to i) poor description of the metal metal bonds and ii) 
basis set superposition effects (i.e. the H basis functions mimic the effect 
of the missing 4p functions and thus lower the cluster energy even without 
the H atom electron and associated nuclear charge). Thus, unlike Upton and 
Goddard [14] who omitted the important 4p functions, we have included the 
4p functions for those Fe atoms which are within bonding distance of the 
H atom in any site considered in the cluster. These Fe atoms are denoted 
as primary (P) . To reduce edge effects we have also added the nearest neighbors 
of the primary atoms. These we denote as secondary (S) atoms. Because 
these atoms are more distant from the H adatom, the omission of the 4p functions 
is less serious and we chose to include only the 4s functions on these atoms 
in order tc keep the calculations of more reasonable size. For example, 
the Fe,, cluster involves 330 basis functions with 4p functions on all centers, 

DO 

but only 198 functions if 22 atoms are primary and the remainder are secondary. 


Fig. 2 shows the Fe^g cluster which is a representative cluster for interior 
sites. The primary atoms here consist of eight atoms in an arrangement consist- 
ing of two fused tetrahedra. Each primary atom has the full complement of 
eight nearest neighbors which are included at the secondary level. The loca- 
tions of the H atoms are also indicated in Fig. 2. The Fe^g cluster contains 
i) octahedral sites which have six nearest neighbor Fe atoms and ii) tetra- 
hedral sites which have four nearest neighbors. (Note that these are not 
regular octahedra or tetrahedra since they are compressed in one direction.) 

Moving from left to right along the h atom positions indicated hi’.PPen circles one 
passes alternately through octahedral and tetrahedral sites at separations of d/4 
where d is the lattice constant. Thus the potential should exhibit a periodicity 
in the direction indicated with X = d/2, i.e. all octahedral sites should be 
equivalent and all tetrahedral sites should be equivalent. This periodicity 
is not necessarily present in the clusters due to cluster edge effects and 
we have used periodicity as a criterion for judging the size of cluster edge 
effects. Two other clusters were studied which are related to the Fe^g 
cluster. These were the smaller Fe^g cluster which is related to the Fe^g 
cluster by deleting two primary atoms and associated secondary atoms from one 


end of the cluster and a larger Fe^g cluster which extends the Fe^g cluster 
by addition of an additional tetrahedron of primary atoms and associated 
secondary atoms at one end. 

The BCC metal also exhibits periodicity in a direction perpendicular to 
the path discussed above for the Fe^g cluster. This path is indicated by the 
solid circles in Fig. 2. Also illustrated in Fig. 2, by the triangle, is a 
trigonal site which is found to be a saddle point connecting the two adjacent 
tetrahedral sites. In order to study the T-<-*T diffusion via the trigonal site 
we considered additional clusters which are equivalent in the two periodicity 
directions. The primary cluster here is a 22 atom cluster consisting of four 
fused BCC unit cells centered about the intersection of the two paths given 


above. Here we considered an Fe^j cluster with all atoms primary and an Fe 
cluster which Includes in addition the nearest neighbors as secondary atoms. 


66 


Fig. 3 shows the Fe^g cluster for the surface 1 sites. Here we are con- 
sidering two-fold and four-fold sites. The Fe^g cluster has six surface and 
two second layer primary atoms. The secondary atoms consist of ten in the 
first layer, six in the second layer, and six in the third layer. Note that 
some of these are at second nearest neighbor distances. The Fe^g cluster 
adds to the Fe^g cluster three additional primary atoms (two in the first 
layer and one in the second layer) and associated secondary atoms leading to 
an additional four-fold surface site. For each of these clusters the H atom 
was moved perpendicular to the surface for the two-fold and four-fold sites. 

The ground states of the clusters were obtained by an Auf Bau method. 

The symmetry used in the calculations is C2 V for the interior site clusters 

and C for the surface site clusters. The orbitals are filled using the Auf 
s 

Bau principle until the full complement of electrons is present. The cluster 
configurations for the interior site clusters are as follows: Fe^g has six a^, 

three b2» three bp and two a 9 orbitals doubly occupied and one a^ and one b2 
orbitals singly occupied which leads to a state. The triplet state arises 
here because the singly occupied orbitals belong to an E representation in the 
full cluster symmetry which is higher than C2 y . Fe^ has eight ap four b2» 
four b^, and two a2 orbitals doubly occupied. Fe^g has eleven a^, five b2> 
five bp and three a2 orbitals doubly occupied. Fe^g has twelve ap nine b2> 
seven bp and five a2 doubly occupied. The cluster configurations for the surface 
site clusters are as follows: Fe^g has nine a' and six a" orbitals doubly 

occupied. Fe^g has twelve a' and seven a "orbitals doubly occupied and one a" 
orbital singly occupied. The Fe^H clusters are all doublet states with the 
exception of the Fe^gH interior sites cluster which is a quartet state and 
Fe^g which is a singlet state. 


7 


IV. Discussion 


Fig. 4 shows the binding energies for the Fe^gH, Fe^gH, anc * Fe^gH c l usters * 
As discussed in Section III, for the solid all tetrahedral sites are equivalent 
and all octahedral sites are equivalent. This behavior is most closely approx- 
imated by the Fe^gH cluster, where we see pseudo periodic behavior with the 
tetrahedral sites slightly lower in energy than the octahedral sites. Note 
that these calculations do not include enough points to demonstrate the shape 
of the potential (i.e. one would expect an approximately sinusoidal variation 
between the calculated points). The other clusters, on the other hand, at 
first appear to exhibit rather different potentials with the Fe^gH cluster 
showing a decrease in energy in moving from the center to the edge of the 
cluster while the Fe^gH cluster shows an increase in energy for motion along 
the same path. However, the potentials can be made to look very similar for 
the three clusters by including a correction term (for edge effects) which 
is linear with distance. Selecting such a correction term to make the two 
octahedral sites degenerate leads to the corrected curves in Fig. 4b which 
show a remarkable similarity to each other. From Fig 3b one sees that in all 
cases the tetrahedral sites are below the octahedral sites. The calculated 
separations are 0.21eV for Fe^gS, 0.30 and O.lOeV for Fe^gH, and 0.28 and 
O.lleV for Fe^gH. Thus, the calculations can be interpreted to indicate 
a tetrahedral to octahedral separation of 0.10 - 0.30eV. 

Fig. 5a shows two cuts through the potential surface for the FeggH cluster. 
One path follows one of the periodicity directions in Fig. 1, while the other 
path is from the central octahedral site toward the trigonal site. Fig. 5b 
differs from Fig. 5a in that a linear correction has been added to make the 
potential flat from the center of the cluster to the edge (as was done for 
Fig. 4). These two cuts through the surface show that the octahedral site 
is a maximum on the potential energy surface and are consistent with the tetra- 
hedral site being a minimum on the surface with the trigonal site being 
a saddle point connecting two adjacent tetrahedral sites. A remarkable 
feature of Fig. 4 is the clearly sinusoidal variation of the energy along the 
periodicity direction. This result suggests that the FeggH cluster is exhibit- 
ing true periodicity although from comparison of Fig. 4a and Fig. 4b there 
still appears to be a significant edge effect. 


From Fig, 5b the tetrahedral to octahedral separation is 0.23eV and 
the tetrahedral to trigonal separation is 0<17eV (based on the central octa- 
hedral site and the adjacent tetrahedral site) . These orderings may be under- 
stood based on the following simple geometric argument. The octahedral site 
does not correspond to a regular octahedron but has two shorter FeH distances 
(2.703a Q ) and four longer FeH distances (3.823a Q ). The tetrahedral site has 
all FeH distances the same (3.022a Q ) and the trigonal site has two short 
FeH distances (2.819a Q ) and one longer FeH distance (3.023a Q ). With the 
one-electron ECP and the same basis set used in the calculations diatomic FeH 
has a bond length of 2.96a Q . Clearly one expects an H atom in the solid to 
prefer an FeH distance longer than for the diatomic molecule. Thus, it should 
be downhill for the octahedral site with two short bonds to distort toward 
the tetrahedral or trigonal sites. The tetrahedral site is favored over the 
trigonal site since it has four reasonable FeH bond lengths while the trigonal 
site still has one compressed bond distance. 

For the octahedral, tetrahedral and trigonal sites some relaxation of 
the lattice was allowed both with and without the H atom present. The clusters 
here consisted of six, four, and three primary atoms plus the nearest neighbors 
as secondary atoms. The resulting cluster sizes are Fe^g, anc * ^ e 20’ re “ 

spectively. Fig. 6 shows the distances which were varied in the geometry 
optimization. Here only the locations of the primary atoms were varied while 
the secondary atom locations were fixed. For the bare clusters the Fe^g 
and Fe^ cluster geometries were optimized with the constraint that r^ and r 2 were 
varied simultaneously; this lead to geometric expansions compared 'to the BCC 
geometry of 3.4% and 3.9% for octahedral and tetrahedral, respectively. For 
the trigonal cluster r^ and r 2 were varied separately giving a 2.4% increase 
in r 2 and a 6.5% increase in r^. For the clusters with the H adatom the in- 
creases compared to the bare clusters were: r^=+10.2% and for octa- 

hedral, rj=+3.8% and r 2 =?+12.2% for tetrahedral and rj=+9.7% and r 2 =+1.2% for 
trigonal. The corresponding increases in binding energy are 0.27eV, 0.22eV 
and 0*29eV for octahedral, tetrahedral, and trigonal. Anplying these correc- 
tions to the energies from Fig. 5b leads to an estimated tetrahedral to 
octahedral separation of 0.18eV and an estimated tetrahedral to trigonal 
separation of O.lOeV. The volume expansion for the tetrahedral site, which 
is the only site which represents a minimum on the surface, is +21%. 
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Fig. 1 shows the potential for moving a H atom perpendicular to the (100) 
surface for two-fold and four-fold sites of the Fe^K surface cluster. Here 
we see that surface atoms are much more stable ('VL.3eV at the four-fold site) 
and have a minimum at 'v0,5a o above the surface., We expect a barrier for 
pushing the H atom beneath the surface for a four-fold site because it would 
be moving in the direction of a second layer Fe atom. For the two-fold sites, 
on the other hand, a location in the surface is basically an octahedral loca- 
tion and from the studies of interior sites we expect it to be downhill to move 
to the adjacent tetrahedral site. Indeed it is found that the H atom will 
move from this location to the interior. Fig. 5 shows similar calculations 
for the Fe^H surface cluster. Here we see that the Fe^H cluster shows results 
which are consistent with the Fe^R cluster. The binding energies here are ^1.7eV 
for site A and njl.3eV for site B. This difference between the two four-fold 
sites is another manifestation of cluster edge effects. 

As discussed above we expect a large barrier to moving an atom from a 
four-fold site directly toward a second layer atom. Moving from a four-fold 
to a two-fold site is uphill by M..25eV but it is downhill from here to an 
interior site; this implies a barrier of 'U).25eV in addition to the endother- 
micity of 'vl.OeV to the process of moving a H atom from a four-fold surface 
site to an interior site via a two-fold site. However, there may be alternative 
pathways such as penetrating the surface at a four-fold site followed by side- 
ways movement toward an adjacent tetrahedral site which involves little or 
no barrier. 

From examination of Figs. 5, 7, and 8 we see that the most stable site 
is the four-fold surface site with an estimated binding energy of 1.3 - 1.7eV 
while the two-fold surface site shows essentially zero binding energy (based 
on the Fe^H cluster) . Interior sites are less stable than surface sites being 
bound by about 1.3eV (based on Fe^H) . This result is consistent with the 
experimental observation that H atoms are more stable at the surface with only 
a small fraction of the H atoms moving inside the solid (see section V) ; 
however, it is not clear that the binding energies are well enough converged 
with respect to cluster size to strongly support this conclusion . 




^ V. Comparison to Experiment 

The current picture of 1^ chemisorption on an iron surface is that the 
molecule first bonds in a weakly bound state which subsequently dissociates 
j. leading to two H atoms chemisorbed on the surf ace [21]. The detailed energetics 

for H atoms interior to the bulk are not known but overall the interior sites 
are believed to be less stable than the surface sites [22], H atoms are known 
to be quite mobile inside the bulk [20], The solubility and diffusivlty of 
hydrogen in iron has been reviewed by Kiuchi and McLellan [24]. The quantities 
which are experimentally accesible (experimental value in parenthesis) are 
i) the heat adsorption (0.9eV) [22] which is the energy for chemisorption of 
on an iron surface (normally as film) . Since is known to dissociate on 
- the surface this quantity is twice the binding energy of an H atom minus the 

D £ of ii) the heat of solution (0.3eV) [23] which is as in (ii) except 

that the two H atoms are in the bulk. Based on an analysis of non-Arrhenius 
behavior of the solubility it has been inferred that tetrahedral sites are 
« about 0.23eV more stable than octahedral sites, iii) the activation energy 

•i 

for diffusion in the bulk iron lattice (O.leV) [20], and t iv) the volume expansion 
for a hydrogen atom in the perfect lattice which is £0.17 atomic volumes for 
BCC iron [25]. 

v Using the experimental of Hg in conjunction with the heat of adsorption 

and the heat of solution leads to binding energies of 2.8eV for one H atom at 
the surface and 2.5eV for one H atom in the bulk. The best calculated values 
are less than this being 1.3-1.7eV (four-fold, surface site) and 1.3eV (inferior 
f tetrahedral site), respectively. This is not an unexpected result given the 

SCF model which typically underestimates binding energies and considering 
that the same ECP and basis set underestimates the diatomic FeH binding energy 
by 1.4eV. 

$ The activation energy for diffusion is estimated to be O.lOeV based on 

the Fe^H cluster. One expects that the effective barrier should be lower 
than this , since tunneling is known to be important for hydrogen [26], Thus, 
it appears that the calculations underestimate this barrier since the experi- 
| mental barrier is thought to be 'V'O.leV. Finally, the calculated volume ex- 

pansion is 212 in reasonable agreement with the experimental estimate of 17%. 
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These encouraging results Indicate that the SCF cluster model In conjunc- 
tion with the one-electron ECP and basis set used here provides physically 
reasonable results and can be relied upon to provide useful information in 
situations which are experimentally inaccesible, e.g. the interaction of H 
atoms with defects in the solid. 
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Table I Parameters of the GTO Fit of the One-Electron 


Fe Effective Core Potential a 
N ZETA COEFFICIENTS 


F Potential 


2 

2 

2 

2 

2 

2 

1 


0.47390366 

1.58405370 

4.95107780 

19.31682400 

86.61302200 

594.26962000 

32.93220700 


-1.26119200 

-9.53136520 

-25.40468800 

-64.39665200 

-190.95396000 

-733.44870000 

- 1.00000000 


S-F Potential 


2 0.77791885 
2 2.35091970 
2 7.40039490 
2 27.76902900 
2 175*84886000 
0 488.84440000 


4,11278120 

20.80044300 

50.55318300 

126.11910000 

515.48734000 

3.00000000 


P-F Potential 


2 0.62727306 
2 1.86402300 
2 5.77181370 
2 18.09789100 
2 66.97594900 
0 143.17037000 


2.35615760 

13.98017900 

50.95290200 

144.04543000 

376.36175000 

5.00000000 


D-F Potential 


2 

2 

2 

2 

2 

0 


0.43787116 

1.90303490 

7.19745380 

24.49263100 

81.70534400 

152.31608000 


1.73163370 

11.35304600 

45.18231400 

172.21949000 

582.84560000 

7.00000000 
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The Fe one-electron ECP was developed by Dr. Louis Kahn [2], 


Figure Captions 


Fig, 1 Calculated potential curves for Fe 0 states arising from the 4s 3d q + 
17 ^ 

4s 3d atomic limit. Curve A is for the high spin state which has no 
3d bonding, while curve B is for the low spin state which could have 
maximal 3d bonding. 

Fig. 2 The Fe^ cluster for interior sites in BCC Iron. The locations of 
the primary atoms are indicated by the large circles; the locations 
of the secondary atoms are indicated by the intermediate size circles. 
The H atom locations for one periodicity direction are indicated by 
the small open circles; the H atom locations for the other periodicity 
direction are indicated by the small solid circles, while a trigonal 
site is indicated by a triangle. 

Fig. 3 The Fe^Q cluster for surface sites. The conventions are as for Fig. 1. 

Fig. 4 Calculated binding energies for octahedral (0) and tetrahedral (T) 

sites in BCC Iron. Fig. 4a shows the calculated points while Fig. 4b 
shows the effect of adding a correction term (see text) to compensate 
for cluster edge effects. 

Fig. 5 Calculated binding energies for two cuts through the potential energy 
surface for H in the Fe^ 6 cluster. Curve A is for motion along one 
of the periodicity directions while curve B is for motion in the direc- 
tion of a trigonal site. Fig. 5b differs from Fig. 5a in that a 
correction term has been added to compensate for cluster edge effects 
(see text). 

Fig, 6 The geometric parameters which were varied In studying lattice relaxa- 
tion for H in octahedral, tetrahedral, and trigonal sites. 

Fig. 7 Calculated binding energies for two-fold and four-fold sites on the 
(100) surface of BCC Iron from the Fe^gH cluster. 

Fig. 8 Calculated binding energies for two-fold and four-fold sites on the 
(100) surface of BCC Iron from the Fe^H cluster. 
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APPENDIX B 


The following papers describe the calculations for the TM 
hydrides and dimers. 


CASSCF/CI calculations for first row transition metal 
hydrides: The TiH( 4 <f>), VH( 5 J), CrH( 6 i:+), MnH( 7 i: + ), 
FeH ( 4 6 J ), and NiH ( 2 J ) states 

Stephen P, Watch*’ and Charles W. Bauschlicher, Jr. b) 

Polyatomic s Research Institute Mountain Heu', California 94043 
(Received 2J October 1982. accepted 28 December J982) 

By consideration of atomic coupling arguments and interaction of 4rM' and is '!</"*• derived terms, the 
ground states of the transition metal hydrides are predicted to be TiH ( 4 4t>), VH |'d). CrH (\F*), MnH (\F*|, 

FeH (‘A I, CoH i‘0l, and NiH | : d ). All of these systems have been studied by a CASSCF/CffSDJ procedure 
with the exception of CoH. The calculated parameters derived from the CASSCF/CI|SDJ potential curves are 
in good agreement with the experimental values where known (no information exists for TiH and VH). 
Inclusion of atomic correlation leads to significantly belter agreement with experiment particularly for R,. 

These improvements are related to a more balanced description of the atomic states in the correlated wave 
functions. 


ORIGINAL PA&E IS 
OF POOR QUALITY 


I. INTRODUCTION 

The chemistry of the transition metals (TM) Is es- 
pecially diverse because of the presence of several 
low-lying atomic states which may be utilized In bonding. 
For the first transition row the lowest atomic states 
are derived from the 4s 2 3 tf , 4s 1 3d 7 * 1 , and 4s l 4/> l 3d" 
configurations. As we will see the bonding in TM In- 
volves a complex Interaction among terms arising out 
of the various atomic levels. A number of studies 1 ” 5 
have shown that the separations of the TM atomic states 
are not well described at the Hartree-Fock (HF) level; 
however, inclusion of electron correlation does sig- 
nificantly Improve the description of the atomic states. 
Thus, electron correlation has to be included in the 
molecular calculations. 

There have been a number of previous theoretical 
studies of the TM hydrides. Earlier all -electron studies 
by Scott and Richards 1 ” 6 and Bagus and Schaefer 8 used 
HF wave functions which completely neglect electron 
correlation. More recently Henderson, Das, and Wahl 
(HDW) 10 and Das 11 have carried out limited MCSCF 
studies of the TM hydrides. In these studies atomic 
correlation was largely neglected. This was justified 
in the work of Das 11 by the observation that for VH the 
calculated properties of several states were similar In 
the GVB/POL-CI calculations of Walch 12 which included 
limited correlation (4s 2 -4/> 2 ), and In the calculations of 
HDW which neglected these correlation terms. How- 
ever, Bauschlicher and Walch 15 pointed out that for ScH 
neglect of 4s*-4/>* near degeneracy terms leads to an 
incorrect ordering of the l E* and 5 A states. Similarly 
for NiH Walch and Bauschlicher 14 found that a correct 
description of the state of NiH required a balanced 
description of the 4s* 3d* and 4s 1 3d* states of Nl, and 
consequently an accurate description of electron cor- 
relation was necessary, This conclusion is also sub- 
stantiated for NiH by the CASSCF/CI calculations of 
Blomberg el al. is and had been suggested by the earlier 

^Srpported under NASA grant NCC2-148. 

^Present address: NASA Ames Research Center Moffett Field, 
California 94035, 

* Mailing addresss 1101 San Antonio Road, Suite 420 Mountain 
View, CA 94043. 


studies of Bagus and Bjorkman 14 and R, Martin 17 In 
which limited 3d shell correlation had been included. 

For completeness we also mention the work of Goddard 
et al, lt on NIH In which electron correlation effects for 
the Nl atom were Included In a semlempirical way using 
a modified effective core potential, This approach led 
to properties for NiH In good agreement with experi- 
ment at the two configuration MCSCF (GVB) level, 

In the present study calculations were carried out for 
the predicted ground states of TiH( 4 $ ), VH( 5 A), CrH( 6 2*), 
MnH( 7 S*), and NIH( 2 A). For FeH both the °A and 4 A 
states were studied, since both are likely candidates 
for the ground state. As discussed In Sec. H, the 
ground state symmetries are predicted based on a com- 
bination of atomic coupling arguments and coupling of 
4s*3<f amd 4s 1 3d n * 1 terms in the molecular system, 
Electron correlation Is included by a CASSCF/CI(SD) 
treatment. The CASSCF Includes near-degeneracy 
effects (4s 2 - 4 {?) while correlation of the 3 d electrons 
Is Included at the Cl level. 

II. QUALITATIVE FEATURES OF THE BONDING IN 
TM HYDRIDES 

Figure 1 shows the relative ordering of the 4s 2 3d" 
and 4s 1 3d"* 1 states for the first transition row. 18 The 
trends in Fig. 1 may be understood in terms of two 



FIG. 1. If* 1 — 4s ! 3d" excitation energies of scandium to 

copper Iftttf* 1 ) -Els 7 d*)]. All units are in eV. 
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TAB LX I. Calculated parameter* (or the TM hydride*. 


91 at e 


R, (A) 

(cm M ) 

0, (#VI 

3 d population 

TlH '♦ 

PW* 

l.M 

1407 

2.12 

2.26 


SR* 

1.94 

1410 




D • 

1.91 

1331 

1.7 


VH ‘X 

PW 

1.74 

1590 

2.30 

3.42 


D 

1.74 

1609 

1.77 


CrH *Z* 

PW 

1.70 

1405 

2.10 

4.0b 


O 

1.71 

1570 

2.0 



Expt.' 

1.66 

1591 



MnH V 

PW 

1.77 

1039 

1.71 

5.06 


O 

1.94 

1432 

1.9 



txpt.' 

1.73 

154* 



F*H ‘a 

PW 

1.72 

1500 

1.95 

6.00 


r 

1.02 

I960 

1.7 



D 

1.77 

1300 

1.43 



SR 

1.73 

1050 




Expt.' 

1.77 .0.05 




FeH *A 

PW 

1.90 

1710 

(3.20) 

6.26 


J 

1.50 

1730 

(2.0) 



Expt.* 

1.03 

1000 



NIH 

PW 

1.47 

1902 

2.79 

0.66 


BSR* 

1.47 

1911 

2.76 



BB' 

1. 40 





M» 

1.50 

1990 

2.34 



O 

1.55 

1917 

2.5 



GWRl* 

1. 45 

1911 

2.76 



Expt.' 

1.41 

1926.0 

*3 10' 



'Preaent work (PW). 

Scott and Richards (SR) Refs. 6-*. 
*D«a <D> Ref. 11. 

•Reference 2". 


u. r: i. r ; 

OF POOR QUALITY 


•Jaffa (J) Ref. J2. 

'Stevens, Llneberger, and Fetgerle Ref. 30. 

•Davis Ref. 31. 

Sllomberg, Slegbahn, and Roos (BSR) Ref. IS Cl results. 0, referenced to 4s 1 3d*. 

'Bag os and BJorkman (BBI Ref. 10. The R, value ta from Cl calculation*. 

'Martin (M) Ref. 17 Cl values. 

'Goddard, Walch, Rappe, and I'pton 1GWRL') Ref. IS. 

l Tbe calculated 0 # values are with respect to 4s 1 3d* while the experimental 0„ la referenced to 4i J 3^. 
The experimental aeparation Is 0.03 eV with 4s 7 3d*. Thus, the appropriate D, value for comparts no 
to the calculation is 3.16 eV (3. 10 eV with respect to 4* 7 3d*i. 


1 


1 


competing effects (i) a general stabilization of 3 d with 
respect to 4s with increasing Z, and (ti) a preference 
for a maximum number of high spin coupled 3 d orbitals. 
For les6 than half -filled shells both effects tend to 
stabilize 4s , 3d"*‘ with respect to 4s*3d" for Increasing 
Z, thus the 4s 1 3d"' 1 -4s* 3d" separation decreases 
mono tonic ally from Sc-Cr. For Mn 4s*3d‘ is favored 
over 4s 1 3d* because of the extra high spin coupled 3d 
orbital. Thus, at Mn there is a discontinuity with 
4s* 3d’ again becoming lower. For Mn-Cu the 
4s 1 3d"*’-4s l 3d’ separation once again decreases mono- 
tonically due to stabilization of 3d with respect to 4s for 
increasing Z. 

In addition to the 4$ l 3d* and 4s 1 3d* 1 atomic configu- 
rations the 4s , 4/> > 3d’ configuration is also important. 
The excitation energy 4s 1 3d" - 4s' 4p‘ 3d* increases 
monotonlcally with increasing 4 from -16 000 cm* 1 for 
Sc to -28 000 cm* 1 for Cu. 1 * 


The bonding of a H atom to a transition metal atom 
may involve either the 4 s*dd’ or 4s 1 3d*" 1 state. For 
the 4s* 3d" state the bonding involves formation of s/- 
hybrids arising from interaction of the 4s*3<r and 
4s , 4/> , 3d" atomic configurations. This leads to two 
orbitals sz and si which have the qualitative character 
4$ *4pa and 4s -4 po, respectively. The bonding orbital 
(two electrons) has the character of a sigma bond be- 
tween Sc (sz) and H (Is) while the si orbital which is 
singly occupied is hybridized away from the bond pair. 
For the 4s 1 3d”' state the bonding Involves formation 
of a simple Sc (4.0-H (Is) bond. We expect the bonding 
here to be stronger than for the 4s*3d’ state since no 
promotion energy is involved. This picture is supported 
by the population analysis (CASSCF/CISD wave func- 
tions) which shows that MnH which is predominately 
4s* 3d* like has valence s tad p populations of 1.04 and 
0.68 indicating a strong admixture of 4s'4p'3d‘ char- 
acter while CrH which is dominated by 4s 1 3d* has 
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valence » and p population* of 0. 77 and 0. 1 1 , 1. 1 . , vary 
little adn..xure of ip character. 

Considering these two bonding mechanism* in con- 
junction with the variation in the fi'M'-O'M" 1 sepa- 
ration (Fig. 1) one expect* predominate!) 4v*3d* like 
character for element* where 4 * * 3,1 * i* well helow 
4»'3J**' le.g., Sc and Mn), predominately 4s'3d**' 
character for element* where 4*' 3d**' i* veil below 
4.«*3 d* le g., Cr and Cu), and mixed character for the 
other element* with especially strong mixing for 
element* where 4*' 3d**‘ it slightly above 4s 1 3d* (e. g. , 

V and Co). These expectations are borne out by the 3d 
populations in Table I where we see TiH t» predominately 
4«*3d*, VH is a mixture of 4s* 3d* and 4<‘3d 4 , CrH Is 
predominately 4* > 3d > , MnH is predominately 4** 3d*, 

FeH '.a is predominately 4**3d*, and NtH is a mixture 
of 4«*3d* and 4*' 3d*. 

As pointed out elsewhere, '* there is an additional 
complication for ScH where there is a competition be- 
tween bonding to s? and 3 Jo for the 4** 3d' configuration. 
However, we suspect that a “d-bonded" ground state is 
peculiar to ScH for the following reasons (i) As one 
moves from the left to right side of the first transition 
row both the 3d and 4» orbitals contract, but the ratio 
of (r*, v’ w / increases monotonkcally from 2.364- 
3. 239.* Thus bonding to the 4<* pair is favored in- 
creasingly as one moves from Sc to elements on the 
right side of the row. (ii) Formation of the “d-bond" 
results in loss of exchange energy as the number of 3d 
electrons increases This effect should be most sig- 
nificant in the center of the row where the maximum 
number of high spin coupled electrons occurs. 

In the previous discussion we considered the inter- 
action between H (ls)and the 4«.*3d" and 4s ' 3d"*' atomic 
configurations. We now concider how to select the 
lowest 3d orbital occupancies. As discussed in some 
detail elsewhere,* 0 just as one may write a given atomic 
state as a mixture of determinants it is possible to 
express a given determinant as a mixture of atomic 
states. For example, for the Ti atom in the 4s* 3d* 
state the configuration 4'*3d' ,, 3d’'' is pure but the con- 
figuration 4s*3do' 3d* 1 is a mixture of 40*4 *F and 60$ 

*P. These relationships have been worked out for the 
3d', 3d*. 3 d\ 3d 4 , and 3d 1 configurations and are given in 
Table U. Given this information one then expects for 
Ti that the 4**3do’ 3d*' configuration is 0.62 eV above 
4s*3d(' 3d&' (i.e. , 0.6« 1.03 where 1.03 is the excita- 
tion energy *F — • *P). Thus In the absence of other 
effects one expects it to be more favorable to have the 
configuration 3da*3db x than 3do' 3d>' and similarly for 
the other TM hydrides one may pick likely candidates 
for the ground states from Table 1. Note that here we 
are making use only of atomic information. 

Now consider forming a TM hydride. As a first case 
consider the stale of ScH arising out of 4<*3d'. The 
SCF configuration here is 

6* s/‘3d' o3oo , (1) 

where b is the bond orbital [sz -HJ si like' and si is an 
oroital of 4» - is character. Looking at Table II we see 
that the three components of *£> are degenerate and no 
direct information as to the favorability of 3 do, 3d*, 


TABLI D. Composition of determinants 
In term* of pure atomic slates. 



3d* 


(21 

100$ 

»0 

(II 

100$ 

*0 

<01 

100$ 

'0 


Id* 


(211 

100$ 

'F 

(201 

100$ 

'F 

(III 

00$ 

*F. 201 l P 

<2ll 

•0$ 

*F . 40$ 'P 

(101 

40$ 

} F. 60$ l P 

(231 

20$ 

*F, 00$ 'P 


3 ^ 


(2101 

100$ 

'F 

(2ll> 

100$ 

*F 

(203) 

•0$ 

*F . 20$ *P 

(2 Oil 

• 0$ 

*F. 40$ ‘P 

(213) 

40$ 

*F. eo< *p 

(toll 

201 

*F. 00$ *P 


3d* 


«ioi» 

100$ 

'0 

U102> 

100$ 

■D 

tfiTTi 

100$ 

»0 

(2oT3i 

100$ 

*0 

(ioT5» 

100$ 

‘0 


a dr 


( 2 loHl 

100$ 

*S 


or 3dt is obtained from purely atomic information. 
However in the molecular symmetry 3 Jo and 4s mix 
and in Eq. (1) allowing si to mix in 3d<r character leads 
to mixing in a piece of 4s' 3d* like character: 

6 particle 6*sl* 3d6* adaa *x b t 3do , 3db x adaa 100>*f; 
vparticle 6*si' 3d*‘ o^oo ♦* 6*3dr J 3dr' opoa 40$ *f 1 
o particle 6* si x ida 1 afiaa . (2) 

Note here that the 3dfl particle leads to 4s*3d<r 1 3db x 
which Is pure l F, the 3d* particle leads to 4s*3do , 3d*' 
which is 40$ *F and 60^ *P, while the 3dr particle leads 
to no 4s '3d* term due to the Pauli principle. Thus 
hybridizing 4s and 3 do Is equivalent to mixing 4s* 3d 1 
and 4s' 3d* character and one expects this process to 
be more favorable for a 3 df particle than for a 3d* 
particle which is in turn more favorable than for a 3 do 
particle. On this basis one predicts the molecular 
ordering *A<*I1 <*£’ which is In fact the calculated 
ordering. “ 

Considering now TiH, from Table n the lowest atomic 
configurations are 


b l si x 3dr‘ 3df' a 

(3) 

6* si 1 3do' 3 db x afioaa 

(4) 


which are pure l F for the atomic case. However, Eq. 
(3) mixes with 6*3do' 3d*‘ 3dt' which Is pure *F, while 
no such mixing is allowed for Eq. (4) Thus one pre- 
dicts the ground state of TiH is 4 4> arising from Eq (3). 
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FIG. 2. Origin of VH valence slalrs. The flrat column show* 
the ordering of the 4s : 3d 3 derived term* and the last column 
show* the ordering of the 4«' 3d* derived term* while the cen- 
ter column them » the calculated ordering ol molecular atatet. 


Similarly for VH the moat favorable choice is 
ft'si'Sdfl'Sdir.'Sdr,' , (5) 

which mixes with 

b*idr l 3ifl' idr 1 3dr ( ' adaoa , (6) 

leading to a 'A ground state. 

For CrH mixing of 4s and 3 do couples 
b l tV 3 Ji n l 3 Jr, 1 3dr,' apaaaao 

and (7) 

<>*3dr' 3 Jb, t l 3dC > ( 2. f i 3dir,' 3dr f ‘ adoooao 
leading to a *E* ground state. 

For a less than half -filled 3 d shell the above analysis 
predicts the ground state spin and spatial symmetry 
independent of whether the molecular state is pre- 
dominately 4** 3d" or 4s 1 3d"*'. On the other hand, for 
a more tfan half-filled 3 d shell the spin multiplicity 
depends on whether the atomic character is 4s* 3d" or 
4s* 3d** 1 like. Here the low spin states are expected to 
involve strong coupling of 4s 1 3d**‘ and 4.<* 3d’ derived 
terms as for the elements with less than half -filled 3 d 
shells, but the high spin states should involve less 
coupling of 4s* 3d" and 4s' 3d** 1 because the high spin 
coupling arising from 4s 1 3d**' is expected to be re- 
pulsive. By an argument analogous to the above (sub- 
stituting holes for particles) one predicts that the ground 
state* of the remaining hydride* are MnHf 1,7 !*!, 

Fe'* ,t Ai, CoH(*'*$), and NiH(*‘'A) where the higher spin 
multiplicity corresponds to 4**3d" and the lower multi- 
plicity corresponds to tc'aj"*'. 

To further confirm the qualitative arguments made 
above we also compare the predicted ordering of the 
low -lying states to the calculated POL-CI ordering" 


for VH. Figure 2 shows the cidec. ar terms arising 
out of the V 4* : 37* and 4* 1 3d* atoau; states (based on 
Table 11 and the experimental atomic separations.") Tn* 
left- and right -hand columns of Fig. ? show the atomic 
4** 3d* and 4«'3d* states while the center column 
(labeled “mixed") shows the calculated molecular 
ordering Here the *6 mixed energy if fixrd at the 
energy of the 4** 3d’ (*4 l term since nc mixing of atomic 
slates is expected for *4 (because only one *6 term 
arises). 

The most rctireable feature of Fig. 2 is a large sta- 
bilization of the *A and 'll states as compared to the 
atomic 4»*3d' separation. For the *A stale this eff. ct 
is due to admixture of *A(2lI) terms as above. (Note 
that we use an abbreviated notation for the 3d orbital 
occupancy here and in Fig. 2, e.f., 'a 2ll> Is equiva- 
lent to Eq. (5> and *4(210!) is equivalent to Eq. (6).] 
This admixture of 4s > 3d* and 4s' 3d' character is 
associated with a decrease of 0. 15e c in H, as compared 
to *4, and is also consistent with the CASSCF Cl 
population of 3.42 3d electrons (Table l>. For 'll the 
situation is more complex in that there are two majur 
configurational mixing effects. The dominant con- 
figuration here is *ni20l), howeier. there if significant 
admixture of *I1(2l2)and '11(2102'. The latter configura- 
tion coupling mixes 4s 1 3d' and 4s' id* character and is 
probably responsible for the decrease ir H, by 0.09u e 
as compared to ‘4. Finally, the *1* state shows signifi- 
cant coupling of *E*(10f ) and *(20?» leading to stabili- 
zation with respect to '4 but no admixture of 4**37’ 
and 4s. ' 3d* character which is consistent with *4 and 
'I* having nearly identical R, values. 


III. CALCULATIONAL DETAILS 

The TM basis sets start with Uachiers (14* V; 57' 
primitive set* 1 augmented with Wachter* two 4 p func- 
tions, the diffuse d function of Hay, ** and a set of f 
polarization functions. The ’"exponents used were Ti 
(1.2). V (1 . 4 ). Cr(1.6), Mn(l.l). Fetf.O". Col2.2>. and 
Ni(2.4). These exponents were selected by linear 
extrapolation of the optimum CI(SD) values for Fe and 
Ni. '•*•* These values are near the optimum values for 
a single and double Cl when the 3« and are cor- 
related.' Since we did not correlate the transition 
metal 3* and 3 /> electrons a better choice of exponents 
may have been those optimized for onlv 3d and 4- cor- 
relation.' However, the 4’ functions are found to be 
unimportant for the TM hydrides'* and the 3d correla- 
tion energy and atomic splitting are somewhat insensi- 
tive to the choice of f exponent. Thus, this defect ir. 
the f basis is of no serious conseauence. This basis 
•et was contracted to (5s ip Idlf' using the general con- 
traction feature of IlGCMOU.” 

The H basis set starts with the (6* 1 primitive set of 
van Duijneveldt*' with two additional diffuse * funcu \* 
added in an even tempered manner, while ih* p func- 
tions were selected as a (211 1 contraction of a four 
term CTO fit** to an STO 2p with an exponent of 1.0. 

(8s ip) (5s 3/)]. The Innermost contracted p function is 
similar in spatial extent to the usual single p function, 
while the outer two p functions are diffuse functions 
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4.* 


•4 

to 

>v 

Configuration 

. ; 16 3* 

V 1 i* * • 

4*. 

*% 

4*. 

C orfficlent 

9 

1 

f 

1 

1 

1 

1 

0 

1 

0 

0. to# 

0 

1 

• 

1 

1 

1 

1 

0 

1 

0 

0. 092 

(• 

I 

M 

1 

1 

1 

1 

2 

1 

0 

0. 0*5 

0 

1 

I- 

1 

1 

1 

1 

0 

1 

3 

0. 0*5 


appropriate to H*. Thi» basis set obtain* all but **0.02 
•V of the electron affinity tE. A. ) of H’. 

In the CASSCF calculation* the Ar core* of the TM 
atom* are kept doubly occupied in all configuration*, 
the active *pace con*i*t» of the orbital* derived from 
the transition metal 4«.4/i, and 3/7 orbital* and the Hil.«) 
orbital. I.#., 6 <t-9c. 3*-4-. and 16. In the above the 
9”, 3-. and 1* orbital* are transition metal 3 ./ like, 6" 
i» the bond pair. 7c i* the *7 orbiul, and 6" Introduce* 
left-right correlation of the bond pair, while 4r intro- 
duce* angular correlation of the bond pair. Note that 
at H *• 6" it a transition metal 4« orbital. 7 o is the 
H 1-1 orbital, and 8c and 4* are transition metal 4/> 
orb.tal* (near degeneracy effect), while a* H decreases 
the orbitals evolve much a* in a 2-1 chemical re- 
action.* 1 Note also that for cases where the 3 d n orbital 
i» not occup.ed <7 » H and VH' the 9? orbital is mostly 
Hi2f *) like and there i* no occupied 3i/^ like orbital. 

The 4r orbital it mottiv H 2f r' like, a result which is 
consistent with the nature of the 9c orbital for TiH arid 
VH. Note that the T.V nvdndes contain a significant 
component of TM" H* character (H population 0.22 for 
MnH CASSCF C1SD wave function). 

Because the CASSCF procedure introduces some 
ambiguities into the calculations (e.g., the 9c orbital 
Is hydrogen like for TiH but is transition metal 3do 
like for CrH), it is difficult to obtain a consistent 
description of the bonding In the TM hydrides at the 
CASSCF level. Therefore, the CASSCF calculations 
were followed by CKSD calculations. The starting set 
of reference configurations for these calculations were 
the CVB configurations 


e^7c* 

ec'Bc^c’J 3d" TM(4s*3d‘) -H , 
8 c* 7c 1 

6c* \ 

ec'Bo 1 j 3d * 1 TMi4 < 1 3d"* 1 ) 4 H 

8c* ' 


18a) 


(8b) 


Table ID shows the configurations in the CASSCF for 
the state of MnH. a case which is predominately 
4s* 3d" like. Here one sees that the dominant cor- 
relation terms are the CVB configurations (8a). For 
the *1* state of CrH wtich is derived from 4« , 3d‘ the 
dominant terms in the CASSCF are the GVB configura- 
tions (8b). For this reason configurations (8a) were 
used as reference conf.gu rat ions for Tiri and MnH 


which are predominately 4 «■ * 3 /* like and configurations 
(8b) were used as reference configurations for CrH 
which l* predominately 4v'3d*' like. However, for VH 
and NiH and for the *1 stale of FeH there is a string 
admixture of 4.«*3d’ and 4* 1 3d ** 1 character, and mi re 
extensive reference lists were used as indicated in 
Sec . IV. 

The CASSCF Cl calculations were carried oul with 
Bigcmoi i*’ -**i l*l u*ing the NASA Ames CRAY IS. 
The calculated H,, D, , and values were obtained via 
a Dunham analysis of the points near H,. In those fit* 
in general only three terms are used. Based on a com- 
parison of three and four term fit* for NiH we estimate 
that use of this small number of computed point* lead* 
to error* of ,0 01a c for fi,, (50 cm' 1 for and 
<0.01 eV for D,. 

IV. DISCUSSION 
A TiH 

Table IV present* 1*1 CASSCF and Cl energies for 
the *4 state of TiH, wmle Table I give* the derived 
potential curve parameters. The reference configura- 
tions for the Cl are configurations (8a). Note from 
Table I that the 3d population is 2.26 which is consistent 
with using a 4v*3d* reference. No experimental infor- 
mation exist* for TiH. however the calculated result* 
are compared to the calculated results of Das (D) ,! and 
Scott and Richards (SR).* - ' Here one sees that the 
present results show a shorter bond length and larger 
binding energy— a result consistent with the higher level 
of electron correlation Included in the present calcu- 
lation. 

B. VH 

Table V ») v e6 the CASSCF and Cl energies for the 
state of VH obtained using Eq. (8a) as reference con- 
figurations. Because of the large amount of 4s 1 id* 
character in the Cl wave function (see Table I and 
VIA) more extensive Cl calculations were carried out 


TABLE IN'. Calculated energies for the *4 slate of TIH. 


H . a t 

c Alter 

Cl 

20.0 

- fc 4‘. 82*. 2) 

-$4‘.93‘9 fc 

4.00 

-•4*. 9*3 19 


3.75 

- *4*. 1*92 04 

- Mf>. 01309 

3.50 

-M‘.993‘5 

— M*. (»16 TO 

3.25 

- * 4 *. 990*,1 

- 849.01521 


Vo 7P 7 


4 09'. 
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TABLE V. Calculated energies for (lie »utr 
of VH. 


TABLE Via. 

•t»u oi cm. 


Calculated energies for the 


n‘ C: 

... h c> 


A. •» 

CASSCF 

Cl 

20.0 

- M3. 404 3* 

-M3. 42* *6 

3.60 

- *43.4*>»*1 

-M3. 310 13 

3.23 

-90.46*72 

-*43.31216 

3.00 

-•43.463 13 

-M3. 30* 3S 


A. •» 

CVB 

Cl 

20.0 

- 1043 *33 14 

- 1043. *0o 02 

3.30 

- 1043. M3 ’3 

-1043. ***‘4 

3.23 

- 1043. M4 34 


3.00 

- 1043. *7* M 

- 1043. *41 3o 


TAB) E VIA. VH extended Cl reaulta. 


Case 

Energy 

34 population 

^ reference 

4 4c, 3 reference 

-*43.31216 

3.36 

•2 

44 c, 5 reference 

-*43.31301 

3. 37 

•2 

Mo, 3 reference 

-943.313 33 

3.39 

•3 

Mo, * reference 

-M3. 314 34 

3.42 

*4 


TABLE VI B. VH Mo. 8 reference 
Cl energies. 


A. a, 

£ 

3. SO 

-M3, sms 

3.23 

-*43.31434 

3.00 

-M3. 31062 


(Table VI). In Table VI the designation 4 An denotes that 
there are four active sigma orbitals 6 ct-9<t. The 5Aa 
designation denotes a more extensive MCSCF in which 
a 10? orbital was added to the active space. The lGer 
orbital is mostly 3 da like and is important for describing 
orbital readjustment effects concomitant with the 
mixing of Eqs. (5) an: (6) In the Cl. The reference 
configurations for the extended Cl calculations are given 
in Table VII. The 3 reference list consists of con- 
figurations (8a) while the 5 and B reference lists incl'-de 
the remaining configurations in Table VC. The 5 and 
8 reference lists were selected based on Cl calculations 
using the 3 reference list with the <Ao and $Ac CAS SC F 
orbitals, respectively. Note that the 5 reference list 
includes additional correlation of the bond pair, mainly 
lingular correlation involving H(2/>), which is no in- 
cluded in the 3 reference list. The 8 reference list 
differs from th* 5 reference list »n that single excitations 
from 7 o into the remaining sigma active orbitals are 
included. 


From Table VIA we see that the extended Cl calcu- 
lations lead to results rather similar to the calcula- 
tions using only Eq. (8a) as a reference set. Calcu- 
lations using a larger V basis set lH< 1 1/* 6./ 2^ »/ 

[8* 7/>4<f 2 f) show a very similar population and thi6 
population is only slightly changed by natural orbital 
iterations. From this we conclude that the strong 
mixing of </* and </* observed lit the current CASSCF Cl 
(SD) calculations is converged with respect to basis set 
and level oi correlation. 

Finally, Tabic VI B gives energies obtained at the 
3A< r, 8 •■efere*!.? Cl level. The potential curve param- 
eter* f t* en i . f.i*> » i ire derived from these energies. 
We net* from W.* I that the potential cur.* param- 
eters obtained here are in gn** wi’h the cal- 

culations oi Das" and the e»n-r r*; 'i \i of V ai> **. 1 " 

C. CrH 

Fn- Crj! ground state .» ;r-»<ely 

4s' ad’ like i UK Ca-'SCI was bared on Eq. f*b' 
(3clrlt active space). Tit CASSCf fallowed by a 

SDC7 using Eq. (8b) as references. [Note that the 
CASSCF in this case has a smaller active space than for 
the other cases. The smaller active space here 
omits angula- correlation effects; however, .nese 
effects are included at the CLSD' level and these cal- 
culations arc consistent with those for the other 
hydrides.] The calculated R, and w, are in reasonable 
agreement with experiment 1 ', the bond length is too 
long by ~0.04 A and a-, is too small by 100 cm 1 '. The 
results of the present calculations (Table I and Table 
VIII) are in good agreement with the results of Das 11 — 
a not unexpected result given that the klCSCF model 
used by Das is the same as ou for this particular 
system, and one expects very Ititle mixing of 4.* l 3<.’\ 

O. MnH 

Table DC gives the calculated CASSCF and Cl energies 
for MnH. The reference configurations for the Cl 


TABLE VII. Reference oooflgv rations tor VH extended Cl calculations. 


Reference 

6c 

7c 

ho 

Bo 

JO 0 

16.V 


3», 

« r . 

3% 

«\ 

/ ( 


1 

0 

0 

0 

0 

1 

1 

0 

1 

0 

1 \ 3 

5° 

1 

2 

0 

0 

0 

1 

1 

0 

1 

0 


' 0 

1 

0 

2 

0 

0 

1 

1 

0 

1 

0 

) 6 1 

0 

I 

0 

0 

0 

0 

1 

1 

2 

1 

0 

\ f 

0 

1 

0 

0 

0 

0 

1 

1 

0 

1 

2 

8 I 

2 

0 

1 

0 

0 

0 

1 

1 

0 

1 

0 

I 

2 

0 

0 

1 

0 

0 

1 

1 

0 

1 

0 

' 

2 

0 

0 

0 

1 

0 

1 

) 

s$ 

1 

0 
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TABU IX. Calculated energies tor the *t* 
stair of MnH. 


X. 

CASSCf 

Cl 

20.0 

- 115< . 373 93 

- 1130.439 99 

3. SO 

- 1130.423 M 

- 1150. 301 13 

3. 25 

- 1199.423 4* 

- 1)50.30 2 47 

3.00 

- 1130. 4l« 3i- 

- 1130.497 0 2 


TABLE X. 
ai EeH. 

Calculated energies for th« ‘a ati 

x, •• 

( ASSC l 

Cl 

20.0 

- 1262.911 01 

- 1263.05. 17 

3.50 

-1991. 999 "3 

- 1263. ’14 20 

3.23 

- 1263.00) 36 

-1263. 11716 

3.00 

- 1262.9W U 

- 1263. 114 30 


are Eq. (8a) which is consistent with the 3 d population 
of 5.06. From Table 1 one sees that the calculated 
R, is *0.04 A longer than experiment while the calcu- 
lated is 90 cm J larger than experiment." Comparing 
our results to Das’," we see that his H, is larger than ours 
by “0.07 A land in poorer agreement with experiment), 
•nd his is smaller than experiment by “110 cm*'. 

G > - 1 n the longer h, and smaller in the calculations 
w Las it is puzzling tr note that his D, is larger than 
our a. Note that in a.1 other cases Das’ binding energies 
are smaller. We suspect that overest imation of D, in 
this case is due to neglect of 4‘*-4/>* near degeneracy 
which is more important for the 4«*3 >i % atomic asymp - 
iCIt than for the molecular region. Another puzzling 
fibiiit hen is Das’ comment that the HF configuration 
is not dominant lor MnH. since we find that in the 
CASbCF the HF configuration is dominant for all R. 

(See Table 111. ) 

E. FeH 

Table X gives the calculated CASSCF and Cl energies 
for the '.a state of FeH. These calculations were 
carried out in the same way as the MnH calculations. 

For FeH more extended calculations were also carried 
out to examine the question of the separation of the *a 
and 4 a states of FeH. The *a state of FeH is dominated 

by 


6o* 7o‘ «e' 1 1 1„* Sr,‘ 3» r ' , (9) 

while we find the 4 a state is a strong mixture of the 
quartet coupling of Eq. (9' and 

6e*9cr 1 ie , I i^i UjSr.'Sr, 1 , (10) 

which is best viewed as the 4 a configuration arising 
from 4s 1 3d 1 . In the CASSCF calculations the 4 a state 
has the 7e and 9«* orbitals mixed, i. e. , 7 o • si - 3 dr 
and 9 ~»sz -3 d~. With these orbitals the 4 a state has 
five dominant configurations which are the reference 
configurations giver, in Table XI. The dominant configu- 
rations for the *a state on the other hand, also given in 
Table XI are more clcarl) defined consisting of Eq. (9) 
plus left -right and angular correlation of the bond pair. 

Table XD gives the calculated energies for the *a 
and 4 a states of FeH using the reference configurations 
given in Table XI. H re we find the *a state is 0.09 eV 
below 4 a at the CliSD) plus Davidson's correction level 
of calculation, while the best experimental estimate 1 * 
places *a “0.25 eV above 4 a. 

CVB *1 -2 calculations with the present basis set 
for Eqs. (9 and (10 at h * 20.0e 0 give a separation of 
1.34 eV as compart J to the experimental l F - atomic 
separation of 0.86 eV. Since the 4 a state is a strong 
mixture of 4*-* 3 /* and 4 *• : 3r/ derived terms, one might 
expect that further improvement in the *F - *D separa- 
tion would depress 4 a with respect to *a. However, the 
same CASSCF Cl lSD> calculation with a larger 
|6« 5; 4 f 1 f) basi c set. 14 which gives a X F - X D atomic 
separation of ! . ;. i eV, leads to the same separation of 
0.08 eV with *a oelow 4 a. and the 3 ; population remains 
the same as with the smaller basis set. This esult 
is consistent with cur stjdies 14 of NiH, where we found 
a correct mixing of 4<*3.*‘ and 4s’ 3./* las reflected In 
the bond shortening ai a level of calculation w here the 
separation was still in error by “0.5 eV. 

From Table XII one sees that there is a differential 
Davidson’s correct. on for the 4 a-*a separation of 
FeH of 0. 18 eV (based on *a at 3. 25a c anc* 4 a at 3.00a 0 ). 
Tl. » large differential Davidson's correction is con- 
sistent with the differential percentage reference in the 
Cl (SD) wave functions (92\ for 4 a and 95^ for *a). 
Because of the size of the differential Davidson's cor- 
rection we suspect that the ground state of FeH is 4 a 
(in agreement with experiment) and that the remaining 
error in our calculated separation is due to a need for 


TABLE XI. iteferenoe configuration;; for FeH extended Cl calculation*. 



bo 

7 0 

So 

9 a 16 

.«'r ,6 t* 

3 r . 

4*. 3 



2 

1 

1 

1 

1 2 

1 

0 1 

l 


0 

1 

A 

1 

1 2 

1 

0 1 

0 

*a 

0 

1 

0 

1 

1 2 

1 

2 1 

0 


0 

1 

C 

1 

1 2 

1 

0 1 

A 

* 


2 

2 

0 

0 

1 2 

1 

0 1 

0 


2 

0 

0 

2 

1 2 

1 

0 1 

0 

‘a 

2 

1 

0 

1 

1 2 

1 

0 1 

0 


0 

2 

*» 

0 

1 2 

1 

0 1 

0 


1 

1 

1 

1 

1 2 

1 

0 1 

0 


°* P °0»‘ 


Q i 


/# 


J C^err Pnys Vol 70 Mo 7 1 Apn *983 


4604 


S P Wilch and C W Bau»chlicher, Jr.: Transition metal hyd'ides 


TAKI Nil. ClktlUtd energies (or the ‘aanJ ‘a atatrt ol Ft H. 






4 A 

A. •. 

CASSCF 

Cl 

CASSv F 

Cl 

JC.O 

-l2t;.lM»l 

• IM.NIIf 4 


- 1263. 003 UJ* 


3. 50 
3.25 

-126.'. 999 75 
- UeS.OOl 3- 

- 1263.114 *3 
(- 1263. 12o69>* 

- 1263.117 97 

- 1262. 980 74 

-1363.10*90 

ORIGINAL H 
OF POOR 0 

1 M I'f 

AUTY 

3,00 

2.75 

- 1262.99* 

<- 12*3. 123 741 
- 1263. 11503 
1-12»3.I2I09> 

- 1262. 979 61 
-1262.969 50 

(- 1263.117 711 
- 1263. 107 69 
<- 1263. 120 3*1 
- 1*63. 102 64 
(- 1263. 116 20) 




‘From Table X usir.4 the C\"B releretue configurations |Eq. (dal). 
•Values in parenthesis Include Davidson's correction. 

'From a GVB ♦! »2 calculauoo tor the 'a slate ol FeH. 


higher excitations in our Cl wave fjnction. 

From Table 1 we see that the calculated R, and *>, 
for the state and the calculated R, for the *.3 state 
are in good agreement w ith experiment. J0,n Comparing 
to other calculated results Das 11 studied only the 
state and finds a longer k, smaller w f , and smaller 
D, than in the present calc ulat.or.s, Jaffe >! carried out 
an ECP Cl study ol 'he *4 ana 4 a states of FeH. His 
results show shorter bor.d lengths for both sutes (a 
result which is believed to be due to defects in the ECP); 
however the other calcinated parameters are tr. rea- 
sonable agreement with the present results. 

F. NiH 

For NiH the CASSCF calculations converged to a 
solution somewhat similar to the 4 a FeH calculations. 
The 3 d population here is “6.2 for R - 2. 75« c . Attempts 
to obtain a solution which was dominated by 4 «• 1 3i/* con- 
verged to the same result even when started from 
4s 1 3d* like GVB vectors. It was found that correlating 
the 3<i<?orbital which stabilizes 4v ‘3d* w ith respect to4s*3d* 
lead to stabilization of me 4s 1 3i* like solution. How - 
ever, for consistency, calculations were carried out 
using the 4s* 3d* like orbitals. The GVB-C1 calcula- 


TABLE X1U. Reference configuration* for MH Cl calculations. 


too 

7 o 

Ac 

9c 

16,2-,; 


3% 

3% 

2 

0 

0 

2 

2 

1 

2 

2 

1 

0 

1 

2 

2 

1 

2 

2 

0 

0 

2 

2 

2 

1 

2 

2 

2 

0 

2 

0 

2 

1 

2 

2 


TABLE XIV. 

Calculated energies for the 'a stale of NiH. 

R.a. 

CASsCF 

CAS Cl* 

SC CP 

3.25 

1507.376 39 

- 1 50". 600 35 

- 1507. 60 fl 47 

3. 00 

1507.377 79 

- 1507. *1021 

- 1507.617 35 

2.75 

1507.373 35 

- 1507.61440 

- 1507. 620 65 

2.50 

1507.359 32 

-1507.610 09 

- 1507.614 40 


*Tbc asymptotic energy is - 1507.517 99 Irom a GVB * 1 * 2 
calculation based on a Arif* * Hllri like GVB solution. 


tions were carried out and tne natural orbitals were 
used in a Cl using the reference cc nf>gurations in Table 
XIII. This calculation leads to a 31 population of “8.7 
indicating strong mixing betw een 4r 3 f* and 4v*3d* at 
the Cl level. Note that the results here confirm the 
conclusion reached in the earlier MCSCF POL-C1 
studies 14 of NiH that MCSCF leads tc rrect mixing 
of 4v*3d* and 4<-’3/*. and that this defect is corrected 
at the Cl level. Correlating the 3 / orbitals .ead» to 
a 4s‘3d* like solution as is the case for HF wave func- 
tions. while the present CASSCF calc a.atior.s tend to 
bias the calculation toward 4- l 3d f . 

From Table I we see that the calcaiated parameters 
for NiH (derived from the energies giver, in Table XlYi 
are in good agreement with experiment. 14 Note that tne 
D, value is obtained by comparison to a G\ B ♦ 1 - 2 cal- 
culation for NiH at Urge R '4« : 3<P like solution). This 
is the most appropriate state to compare it to since the 
NiH wave function near R, is mostly 4« l 3d* like. Note 
also that the calculated R, is slightly smaller than ex- 
periment, a result that r» fleets a bias in the present 
calculations toward 4*> J 3d' like orbitals. By contrast 
the MCSCF POL-CI calculations 14 which converged to 
4s* 3d* like orbitals lead to an R, somewhat longer than 
experiment, a result which reflects a corresponding 
bias toward 4s* 3d* like orbitals 

The present calculations are in good agreement with 
the CASSCF Cl calculations of Blomberg el al ., 14 in 
spite of the omission of 4 f functions in the calculations 
of Ref. 15. This result is consistent with the small 
importance of 4 / functions in ScH. l> By comparison the 
calculations of Das 11 show an R, significantly larger 
than experiment, which we view as indicative of the 
need for additional atomic correlation. Interestingly 
the calculations of GWRL’ 14 which are only two con- 
figuration MCSCF but incorporate atomic correlation in 
a semiempirical way are in good agreement with the 
present calculations. 

V. CONCLUSIONS 

The bonding in the TM hy drides is found to involv e 
strong admixture of terms from the 4‘*3d” and 4 i, 3d" 1 
states of the TM atoms. The bonding in the 4« z 3i/ 
state is found to arise by formation of v/> hybrids «z 
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and s: when ». i* 4 • • 4 .* liw«- and *? is 4* - 4 * like. The 
bond pair Invoiur stngltt pairing th« w and H(l«l or* 
bitals while tin *5 orriu. if siugi occupied .*na hy- 
bridized away Iron: th« i»”>nd pair. The bonding in the 
4»‘ 3Y' 1 state involve* formation of a 4* -Hil« 1 bond 
pair. 

By consideration ol atomic coaflinc arguments and 
interaction of 4 ** 3/ and 4* ! 3 / ”* dern ed lerms. the 
ground states of the TM hydrides are predicted to be 
TiH l 4 4>), VH'ai, CrH* MnHiV). FeH 4 ->. 

CoHl’t), and NiH *^ . All of these systems have been 
studied by a CAflICF Cl <SD procedure with the ex- 
ception of CoH. In addition studies have been earned 
out for the *<a state of FeH. 

We find strong mixing of 4.<*3J" and 4s l 3<i”' for 
VH and NiH where the 4**3i' and 4s‘3J*' 1 states are 
close in energy. TiH. MnH. and FeH f*A) are found to 
be predominately 4‘ l 3.f' like, while CrH is found to 
be predominately 4«‘3f** 1 like. These trends are 
consistent with the ordering of the atomic states. 

The calculated H f and L, values derived from the 
CASSCF Cl (SDi potential curves are in good agreement 
with the available expenmer.tal information (no in- 
formation exists for T.H and VH). Inclusion of atomic 
correlation leads to significantly better agreement 
with experiment particularly for R t . These improve- 
ments are related to a more balanced description of 
the atomic states in the correlated wave functions. 
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ponent 3.75. 

»»A. E. Stevens. W'. C. Lienberger, and C. S. Feigerle (un- 
published results). 

’'S. Davis (private communication). 

S: R. Jaffe (unpublished results). 
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The revultv of G\ b MRSD Cl calculation* for the itate of leH* are combined w ith the retultt of prev»u»CAS SCI 
MKSD Cl calculation* for the 4 A and 4 A itate* of I eH to provide theoretical confirmation of the qualitative argument* 
u*ed bv Steven* I leirerle, and Ltneberyer to assign the loweit itate of I eH a* 4 A The calculated electron afftnitv of I ell 
and the calculated R r and u> ( for the ’a itate of F eH ~ are in good agreement with expenmen tal estimates 


Recently Stevens, Fetgerle and Lineberger (1 ) 
(SFLlhave reps’rted photodetachment experiments 
for FeH' and MnH~ For FeH' the lowest photo- 
detachment transition exhibits a vibrational progres- 
sion while a second transition which is 0 24 eV higher 
exhibits no vibrational progression SFL interpret 
these expenments using a theoretical model in which 
FeH' has the configuration 

b*sf’ 3do* 3dfc ? 3drr* (1) 

where b is a bond pair orbital and s? is an sp hybrid 
directed awa> from the bond pair |2] Removing an 
electron from the above configuration leads to two 
low-lying FeH neutral states (i) ionization of the sf 
electron leads to a 4 A state with the configuration 

b 2 3do 2 3di ? 3dn 2 4 A (2) 

and (ii) ionization of a 3do electron leads to a state 
with the configuration 

b 2 sf' 3do’ 3d5 3 3drr 2 6 A. (3) 

SFL argue that (i) involves ionizing an antibonding 
orbital (if t and this involves a significant reduction 
in equilibnum bond length (/?,) while (ii) involves 
ionizing a non -bonding orbital ( 3do) and should in- 
volve little change in R t . Thus SFL assign the lower 
transition with a long vibrational progression to 4 A 
and the upper vertical transition to 6 A. 

* Supported b> NASA Gram NCC2-I48 

54 


The bond length and vibrational frequence (u^i 
for the 4 A state are known directly from espenment 
(3) SFL were able to estimate R t and uj, for the 6 A 
state of FeH and the *A state of FeH - by companng 
computed and experimental spectra. 

Previously Walch and Bauschlicher ( W B > (2) ear- 
ned out CAS SCF Cl SD calculations for the 6 A and 
4 A states of FeH In the present paper GVB Cl SD 
calculations are reported for the J A state of FeH” . 
The calculated spectroscopic parameters are com- 
pared to expenment in table 1 Here one sees that 
photodetachment to the 6 A state of FeH involves 
little change in R t ( <0.01 A) while photodetach- 
ment to the 4 A state involves a large decrease in R t 
(0 12 A ) The calculated R t and u> e for the 6 A state 
of the neutral and for the 5 A state of FeH " are in 
reasonable agreement with the values obtained by 
SFL. Also the calculated electron affinity (EA) is in 
good agreement with the value measured by SFL 
Thus the present calculations in conjunction w ith 
the previous calculations by WB substantiate the as- 
signments of SFL However the qualitative features 
of the wavefunctions are different from the simple 
picture above in that the 4 A state of FeH and the 5 A 
state of FeH' show strong mixing of 3d 6 and 3d 
(see table 1 ). whereas the analysis of SFL is based on 
pure 3d 7 character for these states. 

The basis set used for FeH * is Vi'achters' ( 1 4s q p>d » 
set [4] augmented with (i) two diffuse 4s functions 
selected in an even-tempered fashion, (ii) Wachters' 

0 009 2614/84 S 03 00 (L Flsevier Science Publishers B V 
(North-Holland Phy sics Publishing Division) 


Volume 103. number 1 
Table 1 

Calculated CAS SCI Cl result* for f eH 

ORIGir u } . 
OF POOR QU. 

( HF MICAl PHYSICS IFTTFRS 
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State 


K t t A l 

w, (cm* 1 ) 

T t leV) 

3d 

population 

1 eH *A 

calc 1 2 1 
evp | 1 | 

1.72 

1.77 * 0.05 

1560 

1550 

-008 
0 25 

6 08 

FeH 

calc |2| 
evp | 3) 

1 60 
1 63 

1710 

1680 

000 

000 

6 26 

1 eH* 

calc ■> 

1.72 

1282 

-0 78 

6 50 


evp |l| 

1 79 s 0 03 

1300 : 140 

-093 



*' Present work The calculated spectroscopic constants are obtained sia a Dunham anal) sis of the Cl p'us Dasidson's correction 
energies at 3 00. 3.25. and 3 50e 0 The Cl enerities at 3 25« 0 are — 1 2fc 3 12366 (-1263.14914) where the value in parentheses 
includes Davidson’s correction 


two 4p functions and one even-tempered 4p diffuse 
function, (tii) the diffuse 3d function of Hav (5) . and 
(iv) a single set of 4f functions (o 4) ■ 2.0). The re- 
sulting basis was contracted using the general contrac- 
tion feature of BIGGMOLI |6) (I6sl2p6dlf) 

(7*5 p3d 1 f j The H basis set was the (8s4p) |5s3p] 
basts used in ref (2) The calculations were earned 
out using BIGGMOLI SWEDEN V 

Before considenng FeH “ , we first consider the 
calculated EA for the Fe atom as a function of basis 
set and level of electron correlation (table 2) Two 
different t> pes of Cl calculations were earned out 
The MRSD Cl has as reference configurations the SCF 
configuration for the 4s 1 3d' state of Fe and the SCF 
plus three components of the 4s 2 -» 4p 2 near de- 
generacy for the 4s 2 3d' configuration of Fe" All 
single and double excitations from the above con- 
figurations are included in the Cl The core -valence 
Cl (CV Cl) |7J has the same set of reference con- 
figurations but allows onl> a single electron outside 
the 3d 1 configuration. The qualitative idea is that 
the major differential correlation between the 4s 2 3d' 
and 4s 1 3d 7 states includes the valence correlation 
effects due to the 4s electrons and the core- valence 
correlation effects between the 4s and 3d electrons. 
The intrapair 3d correlation, on the other hand, re- 
mains relatively constant between the two states since 
they have the same number of 3d electrons This ex- 
pectation is substantiated by the results in table 2 
where one sees that the CV Cl result is in reasonable 

* SWLDLV u a vectorized SCI MC SCI direct Cl w ritter 

b> P L.M Siepbahn, B RoovandCW Hausthluhci Jr 


Table 2 

Calculation of electron affimt) (or Fe 


Baviv set 

Method 

Tlectron affimtv ** 

|7s5p3dlf) 

CVCI 

0 72 


MRSD 

0 65 (0.78) b* 

|7s5p4J2f) 

CVCI 

0 76 


evp. f > 

1 00 


»>4s J 3d 7 , '4s 1 3d 7 . 

\alue in parentheses includes Davidson's correction 
c * I \t updated based on the error for CV Cl for Cu and Vi 


agreement with the MRSD * Davidson's correction 
(8] result The (7s5p3dlf) basis set is the same basis 
set used for FeH while the [7s5p4d2f] basis uses 
a more flexible (4d) contraction and replaces the 
single set of 4f functions b> a (3)/ (2] 4f basis based 
on a three-term fit [9) to a Slater 4f with exponent 
2.25 This exponent was optimized for the 4s 1 3d. 7 
state of the Fe atom at the CV Cl level From table 2 
one sees that at the CV Cl level the smaller basis set 
is within 0.04 eV of the larger basis set result for the 
EA 

No experimental value exists for the EA of Fe at- 
om. The EA reported by Hotonand Lineberger (10] , 
0.25 eV, is based on SCF calculations plus an esti- 
mate of correlation effects. A more reliable estimate 
may be obtained by noting that CV Cl with large ba- 
sis sets (equivalent to the larger basis set used here) 
underestimates the E A of Cu (11] and Ni (12] by 
0.23 and 0.24 eV. respective!) Assuming the same 
error for the E A of Fe leads to 1 .00 eV with respect 
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Table J 


Reference ocnfif uratieni for FeH Cl otlculatlom 



to 

7* 

So 

9o 

>***.>* 

»•*> 
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Jo, 
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2 
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0 

1 
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2 

1 

1 

2 

1 

0 

1 

0 
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0 

1 

I 

2 

1 

2 

1 

0 
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1 

0 

1 

1 

2 

1 

0 

1 

2 

I*H 4 a 
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2 

0 

0 

1 

2 

1 

0 

1 

0 


2 

0 

0 

2 

1 

2 

1 

0 

1 

0 


2 

1 

0 

1 

1 

2 

1 

0 

1 

0 


0 

2 

2 

0 

1 

2 

1 

0 

1 

0 


1 

1 

1 

1 

1 

2 

1 

0 

1 

0 

EeH * 

2 

1 

0 

2 

2 

1 

1 

0 

1 

0 


0 

1 

2 

2 

2 

1 

1 

0 

1 

0 


2 

1 

1 

1 

2 
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1 

0 

1 

0 


2 
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0 

1 

2 

1 

1 

0 

1 

0 


to 4s* 3d 7 or 0.1 2 eV with respect to 4s 2 3d 6 . For the 
H atom the basis set used here leads to an error of 
*0.02 eV in the E.A 

The calculations for FeH" used aGVB wavefunc- 
tion to define the orbitals this calculation was fol- 
lowed by multireference singles and doubles Cl 
(MRSD Cl) Given the sue of the Cl calculations the 
use of GVB orbitals for FeH ' as compared to CAS 
SCF orbitals for FeH |2] is expected to be of little 
importance In the GVB MRSD Cl calculations the 
Ar cores of the transition metal atoms are kept dou- 
bly occupied m all configurations, the active space 
consists of the orbitals denied from the transition 
metal 4s. 4po. and 3d orbitals and the H( Is) orbital, 
i.e. bo -9a. 3rr. and 16. In the above the 9 o, 3rr,and 
16 orbitals are transition metal 3d like. 6o is the bond 
parr, 7 o is the sf orbital and 8o introduces left-nght 
correlation of the bond pair Table 3 gives the ref- 
erence configurations which were used in the MRSD 
Cl for the state of FeH The GVB wavefunction 
1 13] includes the first two configurations in table 3. 
The first configuration for FeH" is the SCF configura- 
tion which corresponds to eq (1 ), while the second 
configuration introduces left-nght correlation of the 
bond pair The remaining two configurations corre- 
spond to 3d 6 components of the wavefunction which 
are found lobe sufficiently important that they are 
included as reference configurations in the MRSD Cl. 

The molecular constants given in table 1 were ob- 
tained via a Dunham analysis of the Cl SD plus David- 


son's correction energies As noted in ref. (2) since 
only three computed points are used in the Dunham 
analysis, there may be small uncertainties introduced 
in the denved spectroscopic parameters These errors 
are estimated to be <0 01 A for R t and 50 cm' 1 for 

w e- 

The computed FeH E A is 0 43 eV These systems 
have large Davidson's corrections 0 69,0 35. and 
0.16 eV for FeH" FeH and FeH 6 .i. respec- 
tively These large Davidson's corrections correlate 
with the 7r reference in the Cl wavefunctions 88. 92, 
and 957 for 5 *i. 4 ^. and 6 ^. respectively. Including 
the differential Davidson's correction leads to 0.78 eV 
as our best estimate of the E.A of FeH This value is 
0.15 eV smaller than experiment This result is some- 
what better than might be expected based on the cal 
culated Fe atom EA.but is consistent with the pop- 
ulation analysis which indicates a large amount of H~ 
character Thus the error in the FeH EA may be in- 
termediate between the errors for the EAs of the Fe 
and H atoms. 

We conclude that the present calculations in con- 
junction with the previous results of WB substantiate 
SFL's assignment of as the lowest state of FeH 
with 6 A 0.25 eV higher. (Calculations [2] still show 
the 6 A state 0 08 eV below 4 4. i.e an error of 0.33 
eV in the 4 *J- 6 JL separation ) The key feature here is 
that photodetachment from the state of FeH' 
to give the state of FeH is calculated to lead to 
little change in R t ( <0.01 A) while photodetachment 
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to the *1 state of FeH leads to • large decrease in R t 
(0 12 A) Thus the calculations are in agreement with 
observation of a vibrational progression for photo- 
detachment to but no vibrational progression for 
photodetachment to as predicted by SFL based 
on a simple theoretical model However the qualita- 
tive description of the s>stem is more complex than 
the simple model adopted b> SFL in that the 4 .3 state 
of FeH and the '.3 state of FeH * show strong mixing 
of 3d 6 and 3d' character whereas the analysis of SFL 
is based on pure 3d' character for these states. 

The calculated R t and cc f values for the 4 .3 state 
of FeH are in good agreement with the values mea- 
sured by Davis |3) while the calculated R t and u> e 
values for the 6 3 state of FeH and the *.3 state of 
FeH" are in reasonable agreement with the values 
obtained by SFL from companson of calculated and 
experimental spectra Finally the calculated EA of 
FeH" is in reasonable agreement with the value re- 
ported by SFL. 
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CASSCF CI(SD) calculation* have been earned out for the *12. 1 ij. J i J and 5 A U ttate* of Sc; uun* larjr rauuian 
ba»i* teu The ’ij. 1 Ij and * i ^ Mate* am* from the *D(4i*M' ) ♦ ‘D(4»' 3d 1 1 limit of Scj and are found to be only 
weakly bound (D t * 0 06 e\ and R t • 8 0«o> The * A u ttate arnet from the J D(4t , 3d l ) ♦ 4 F <4i' 3d * 4p 1 1 atomu limit 
Thtt ttate it found to be month bound relative to n* limits ( D t * 0.8 eV and R t * 7 0«o>- 


I . Introduction 

Recently the Sc- molecule has been the subject of 
considerable theoretical studs [ 1 .2] The interest in 
this system anses because of the relevance of metal 
metal bonds to problems in matenals science and be- 
cause theoretical studies find only weakly bound 
states (at least out of the lowest - D(4s- 3d 1 1 ♦ 
2 D(4» 2 3d* ) atomic limit) for Sc- in contrast to mass 
spectrometric experiments [3] which had b*en inter- 
preted to indicate strong bonding 

Das [ 1 ] has earned out an MC SCF study of the 
singlet and tnplet states ansing out of the - D( 4s* 3d 1 ) 
♦ *D(4s*3d' ) and the *D(4s-3d' > ♦ 4 F(4s' 3d* 1 ) 
asymptotes of Sc 2 using a Slater bass set in conjunc- 
tion with a pseudopotential For the 1 1* state of 
Sc 2 Das finds a D e of * 4 kcal mole and an R t of 
*9.5fl 0 which he attnbutes to van der Waals terms 
arising primarily out of the 4s* -* 4p 2 near-degeneracy 
effect for the 2 D(4s- 3d 1 J state of Sc However, as 
Das points out, by making choices for R t and the de- 
generacy factors consistent with the calculations, the 
mass spectrometnc experiments are consistent with 
a binding energy in the range of 3-5 kcal mole. Das 
also considered the possibility of bound states ansing 
from the *D ♦ 4 F limit but found no significant bind- 
ing at least for singlet and tnplet states (quintet states 
also anse from this limit but were not considered) 

Wood, Doran, HilLer and Guest (WDHGi [2] also 
earned out MC SCF Cl calculations using a small 


gaussian basts set. They concluded that the lowest 
state was 5 1 J ansing from the - D ♦ 4 F asymptote 
Their calculations showed a binding energy of I 12 
eV with respect to *D ♦ 4 F and 0.55 eV with respect 
2 D*-D However, due to basis set deficiencies they 
find a 4 F- 2 D atomic separation of 0.5* eV com- 
pared to an experimental separation of 1 43 eV. Cor- 
recting for the error in the asymptotic separation the 
5 XJ state is unbound by 031 eV with respect to 
2 D ♦ 2 D and it does not appear that •* Ij is the 
ground state of Sc 2 as suggested by WDHG 

The present calculations were undenaktn to deter- 
mine (i) the nature of the bonding and the magnitude 
of the binding energy for states ansing out of the 
2 D(4s 2 3d* ) ♦ 2 D( 4s 2 3d 1 ) asymptote, and (u) the 
nature of the bonding ansing out of the 2 D(4s-3d 1 ) 

♦ 4 F (4s 1 3d 2 or 4s 1 4p* 3d* ) atomic limit. The calcu- 
lations reported here are more extensive than those 
of Das or WDHG and involve extensive MC SCF or 
CASSCF followed by CI(SD) using large gaussian 
basis sets. In agreement with Das we find that the 
states arising out of the 2 D + 2 D limit are weakly 
bound, however we use the interacting correlated 
fragments (ICF) method [4] to more accurately de- 
termine the binding energy. We have also investigated 
a S A U state arising out of the 2 D(4s 2 3d * ) ♦ 
4 F(4s 1 4p , 3d* ) atomic limit Here we find a binding 
energy relative to the limits above of * 0.8 eV at an 
R t of *7.0 a 0 We also estimate a binding energy 
relative to the corresponding limit of >03 eV for 
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the state arising f»om *D ♦ 4 F (4s 1 3d* I For Sc 2 thu 
latter limit it 1 43 eV above -D ♦ J D and thu state 
dors not cross the states anting out of *D«*D How 
ever. we And similar states for V 2 and it appears that in 
the case of\' 2 the lowest bond state is of this character 
(5) and thus this state of Sc 2 is relevant to bonding 
in the lowest states of other transition metal dimers. 


2. Caiculationa) details 

The Sc basis set stans with the Wachters ( 1 4s9p5d ) 
primitive set [6| augmented with Wachters' two addi- 
tional 4p functions, thediffuse d function of Hay (7], 
and a set of f polarization functions (a * 1 .4). Two 
different contracted basis sets wc r e constructed from 
this primitive set Basis set I was contracted in a seg- 
mented fashion for use with MOLECULE (8) The 
contraction scheme used here was contraction 3 of 
table VI of ref [bj This contraction was used for the 
s and p functions, while the d functions were con- 
tracted (3111). In this basis set the f functions were 
omitted and the 4p functions were scaled by (1.5 1 1 2 . 
Basis set I was used for the CASSCF calculations 

Table 1 

Reference configurations for the CltSO) calculation 


The final baas is ( 1 4sl lpbdl [8s6p4dj Baas set II 
was generally contracted for use with BIGGMOLI |9) 
resulting in a( 14sl lpbdIO |'sbp3dlf) contracted 
set Here the 4p functions were left unsealed. Thu 
larger basis set was used in the ICF calculations. 

The CASSCF Cl calculations were earned out with 
MOLECULE (8) -SWEDEN [10) using the NASA 
Ames CRAY IS The ICF calculations were earned 
out with BIGGMOLI (4) noname [II] using the 
NASA AmesCDC 7600. 

In the calculations the 18-electron Ar cores were 
kept doubl> occupied in all configurations In the 
CASSCF calculations these orbitals were optimized 
at the CASSCF level, while in the calculations using 
noname the core orbitals were obtained from a HF 
calculation on the 5 1* state of Sc 2 . 

3. Discussion 

The low -lying states of the Sc atom are the 
-D(4s 2 3d’ ) state which u the ground state, the 
4 F(4s’ 3d 2 ) state which u at 1 43 eV. and the 
4 F(4s‘4p , 3d 1 ) state which is at 1.96 eV (12) (the 
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energy levels are averaged over the values I As in- 
dicated earlier the present calculations concentrate on 
states anaui; out of the *D ♦ *D and : D ♦ 4 F limits 
In these calculations we include the 4s, 4p, and 3d 
electrons in the valence space thus the 4s : -* 4p* 
near degeneracy |13) which is important for the 2 D 
state is included at the MC SCF level Correlation of 
the 3d orbitals on the other hand is treated at the Cl 
level 

Table I shows the dominant configurations for the 
J 2" 1 1* J 2* . and ■ ^ states which were consider- 
ed here. The reference configurations in table I are 
the dominant configurations at all R as derived from 
small Cl calculations These configurations may also 
be understood by considering the degeneracies 
which occur at large R (due to the degeneracy of 
the corresponding g and u orbitals at large/? The 
state requires onlv one reference configuration, 
the *1* and *Xr stares require two reference 
configurations, and the 5 ^ u state requires eight 
configurations to dissociate to - Dt 4s’ 3d 1 > ♦ 

4 F (4s- 3d 1 4p* l The remaining two configura- 
tions describe the state arising from - D(4s-3d* ) 
♦ 4 F(4s>3d*l 

In order to keep the CAS SCF calculations of reas- 
onable size, restrictions were placed on the orbital 
occupancies as indicated in table 2. These restrictions 
arise by noting the number of sigma electrons asso- 
ciated w ith each state and allowing up to double exci- 
tations from the sigma block to the pi block These 
excitations correspond to 4s : -* 4prr 2 atomic excita- 
tions. Note that fer the S A U state an additional re- 
striction that the sigma plus pi blocks contain five 
electrons could have been imposed. The CAS SCF cal- 
culations in each case were followed by CI(SD) calcu- 


lable 2 

Constraints on orbital occupancies for the CAS SCF calcula- 
tions 


State 

Constraint 

Number of 
configurations 

*A U 

of iru block 3 -5 fleer ions 

5430 

•iS 

Bfrru block 4 -6 electrons 
np-Tu * pi block 6 electrons 

651 

,£ » 

sip rru block 2-4 electrons 
ngrru * pi block 6 electrons 

2520 


lations from the reference configurations in table 1 

The calculated ene'gies for the C AS SCF C|<SDl 
calculations are given in table 3, while fig 1 shows 
the calc .dated potential curves |C AS S( F CI(SD) cal- 
culations! lor the 1 X* . 3 1* and * ,\ state. 
Qualitatively the 1 2* and - X* states anse from 
2 D(4s 2 3do* ) ♦ ^Dlds^do 1 ) with singlet and triplet 
pairing of the 3do orbitals, while the J 2 f state arises 
from ; D(4s*3drr l ) * • Dt 4»- 3dn 1 ) and the ' Jk u state 
arises out of the 2 D(4s-3do' ) ♦ 4 F(4s l 4po 1 3d6' I 
limit. 

Concentrating first on the states which arise out of 
the *D ♦ *D limit, the basic feature of the calculated 
curves is very weak binding at large /? (* 8 0a 0 ) and 
repulsive behavior at smaller R The weak binding at 
large R appears to anse from the 4s 2 - 4p 2 near de- 
generacy. i.e 4s t -* 4p k X 4s,-*4p f The binding energy 
at theCI level for/? » 8.0a 0 is *0.045 eV for the 1 2J. 
state, * 0 045 eV for the state, and * 0 029 ev 
for the -2* state Note that the CAS SCF curves are 
repulsive at all R Note also that in the van der Waals 
region (i e. >8.0 J 0 )the states arising from the var- 
ious occupancies of the 3d electrons are very close in 
energy. However, at small R the 3d orbital occupan- 
cies do split the curves significantly At small R one 
sees that 1 X* is significantly below 3 2* indicating 
some 3do-Jdo bonding character However it ap- 
pears that the repulsion between the 4s 2 pairs be- 
comes large before sufficient 3do-3do overlap is ob- 
tained to stabilize this state. Given that we are in the 
small overlap region it appeared likely that a one-elec- 
tron bond which varies with distance like 5 would be 
more favorable than a two-electron bond which varies 
with distance like 5v The 3 2^ state ma> be thought 
of as having two one -electron r bonds and indeed at 
small R this state drops below 1 2* although 1 1* and 
3 Xg are nearly degenerate at R values near the van 
der Waals minimum. 

Looking now at the 5 ^ u state, the bonding here 
arises basically out of the 4s 2 + 4s 1 interaction. Based 
on our studies of the transition metal hydrides [14] 
we expect this to be an attractive interaction. The 
reason for this attractive interaction is the presence of 
the low-lying 4s 1 4p> 3d" states of the transition metal 
atom (eg in Sc 4s'4p ! 3d 1 is at 1 .96 eV’y Interaction 
of 4s*3d" and 4s ! 3d"4p I leads to the formation of 
sp hybrids one of which has the character s * z denot- 
ed by sz and the other of w hich has the character s - z 
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Table J 

Tabulated CASCItSDl energyrt (or *in> ut ttatet of S<; •> 


State 


A • 6 0 

A • 6 5 

A • 7 0 

A • 7.S 

A • 80 

A • 1000 

2»A U 

CAS*> 

cl 

-0 36792 

-0 37342 

-0 37616 

-0.36977 

-0 36927 

-0 36125 

*4 U 

CAS 

a 

-043128 
-0 44667 

-043829 
-0 45304 

-0 43943 
-045 344 

-043727 
-0 4 5071 

-0 43346 
-044652 

*1337 

-042288 


CAS 

Cl 

-0 44923 
-0 46754 


-046732 

-0.48270 


-0.47295 

-048575 

-0 47595 

-0.48466 


CAS«» 

Cl 

-0 47J7I 


-048488 


-048632 

-0 48463 


CAS 

Cl 

-0 46023 
-0 47934 


-047010 
-0 48514 


-0 4 7 346 
-0 48594 



In hartree relative to -1519 0 
b) Hi# her root from the * £ u CAS calculation 
c * *Su v «c*or» were used for the Cl 



Fif 1 . Potential energy curves for the low 4yin# states of Scj 
from CAS SCF Cl (SD) calculations The locations of the 
>D« J D and 2 D ♦ 4 F(4s 1 4p l 3d' ) asymptotes are indicated 

denoted by sz This leads in the case of the transition 
metal hydride to an sz -H( It) bond and a singly occu- 
pied sz orbital, and in the case of the transition metal 
dimers to states which may be represented by the fol- 
low ing structures 


O O ■ O t 0— — (I KD 

In the case ot the S 4 U state the remaining electrons 
are in 3do, 4po. 346^ v orbitals which are orthogonal 
to the bond orbitals and therefore the high spin (i e. 
quintet state) is favored. 

From fig. I one sees that the state is bound 
by *0.8 eV with respect to 2 D(4s 2 3d* ) ♦ 
4 F(4s , 4p , 3d* ).The atomic asymptote here is calcu- 
lated to be at 1 .68 eV ( 1 .96 eV experimental I Note 
that experimentally the 4 F (4s 1 3d* ) state is below 
the 4 F (4s 1 4p 1 3d * ) state, while in the calculations 
the 4 F (4s 1 4p* 3d 1 ) state is lower at both the C AS SCF 
and Cl levels This results from the larger 3d correla- 
tion [ 15] in the 4 F ( 4s ' 3d 2 ) state Since 3d correla- 
tion is not included in the C AS SCF calculations the 
4 F (4s* 4p* 3d 1 ) state is lower at the CAS SCF level 
and the resulting orbital bias makes this state also 
at the Cl level This problem makes it difficult to 
study the corresponding 4s 1 3d 2 denved state. How- 
ever. this state is observed as a higher root in the 
CAS SCF (2 5 A g state in table 3) and using these Cl 
energies the binding energy is predicted to be <0.3 
eV smaller for the 5 A g state derived from 2 D ♦ 

4 F (4s* 3d 2 ) than for the 5 A g state derived from 2 D 
♦ 4 F(4s , 4p l 3d* ). Thus we predict a D f of >0.5 eV 
for the S A state derived from 2 D ■*■ 4 F(4s* 3d 2 ) 
Assuming this and noting that the 4s 1 3d"* 1 -4s 2 3d' 1 
separations for Tr and V are 0 81 and 0.25 eV, res- 
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pewtivcl> 1 1 2 ) we predict that the mixed itate * ill 
be unb and b> <0 3 eV with respect to the 4s*3d'' 

* 4i : 3d" limit for Tij and bound b> >0.25 eV for 
the omr bmit in Vj Indeed the lowest stair in Vj 
is found to ante from tmt atomic limit [5] At the 
other end of the row the mued Mates are expected 
to be important for Fe 2 and Co 2 , 

Companni to WDHG. their calculations show a 
deeper well and shorter R t for the 3 I J state {R t • 

4 860 q and O t • 1.12 eV) than do our calculations 
for the state (A, * 7 .O 0 o and D t * 0.8 eV). VXe 
have also studied states of 5 Ij symmetry at the 
MC SCF POL Cl level (16) and find two states, one 
of which is derived from * D(4s 3 3drr* ) ♦ 

4 F(4s , 4po l 3dw* Land the other of which is derived 
from J D(4s*3d6 1 1 ♦ 4 F < 4s* 4po' 3d6* ) Both states 
are bound by *0.8 eV with the first state having an 
R t of *6 5j 0 and the second state having an R t of 

* 7 .O 0 o It is not clear why the results of WDHG are 
so different from ours. We do note however that their 
total energies for Sc 2 are more than 3.0 hartree more 
positive than ours and suspect that some of the differ- 
ence ma> arise from basis set deficiencies in the calcu- 
lations of WDHG (eg basis set superposition error). 

Now w e consider the calculation of the binding 
energy of the states arising out of the - D(4s* 3d * I * 

* D(4s*3d l i limit In these calculations we used basis 
set II The calculations consisted of MC SCF followed 
by POL Cl (I - ) and ICF [4] calculations The ll F 
method has been shown to lead to reliable binding 
energies for weakly bound systems, e g Be : [4) 

Since the minimum in the van der Waals binding 
curve for Sc ; occurs at R > S.0a 0 and the 3d interac- 
tions are small at this distance, we arbitrarily studied 
the 3 1’ and 1 states which arise from 4s 2 3do* ♦ 
4s*3do' (see table 4). 

In these calculations the orbitals were obtained 
from MC SCF calculations with the 4s. 4p and 3do in 
the valence space The MC SCF configurations con- 
ssted of all angle and double excitations within the 
valence space from the reference configurations in 
table 1 (one for 3 1* and two for 1 1* ) These MC SCF 
calculations were then followed by Cl. The POL Cl 
calculations include single and double excitations 
from the reference configurations but with no more 
than one electron outside the valence space To ap- 
proximate the ICF method the orbitals weie localized 
using the transformation 
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Table 4 

Calculated MC SC I POL C I enerpe» for th* ’l J and 1 
states of S<] 


R 

MCSCI 


POL Cl 


lr» 

•u 

_%* 

*1* 


30 0 

-0 50499 

-0 50499 

-0 50*)) 



•0 501)0 

too 


-0 50479 


-0 509)5 

9 0 


0 50414 


-0 50974 

• 0 

-0.50)90 

0 5044) 

-0 509)4 

-0 5091) 

7 5 


-0 50317 


-0 50954 

70 

-0 49415 

-0 50037 


-0 50*39 

to 


-0 4I77J 


-0 49*94 

•• In hartree relative to -1519 0. 


*s • 2 -* 2 *, * 2-> *o u . o, • - 2- ,/J o u . 

and in addition to the POL Cl list configurations of 
the form 

♦i "* v X 0 f - v . 

which allows two electrons into the virtuals were also 
included. (Note that this calculation was earned out 
for the 3 1* $tate since for this state the energy is in- 
vanant to the transformation given above.) F.om 
table 5 we see that the ICF method gives a binding 
energy of 0 046 eV (1.1 kcal mole) while the POL Cl 
method leads to 0 041 eV (for the [6s6p3dlf] basis 
set). Combining this result with the POL Cl estimate 
of the separation between 3 1* and 1 S’ leads to an 
estimated binding energy of 0.061 eV 1 1 4 kcal mole) 
for the 3 Ij or 1 S* state of Sc 2 . 

Das unng an MC SCF model obtains a D t of * 4 
kcal mole, whereas we obtain no binding at the 


Table S 

Calculated MC SCF Cl enerpei for the *1 J state of Sc; 


R 

MC SCF POL Cl 

MC SCF ICF 


|4s6p3d) 

|4t4p)dlf| 

|4i6p3d] 

(4s6p3dlf ) 

30.0 

-0.50*33 h 
OOOOeV 

-0 50947 
0 000 

-030835 

0.000 

-0 50953 
0 000 

1.0 

-0.509)4 

-0.037 

-0 51099 
-0 04) 

-0.51000 

-0045 

-0 51)20 
-0 046 


•) In hartree relative to - 1 5 1 4 0. 
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C AS S( F level Das MC SCF diffcis from ours in that 
he include! onl> the dispersion terms 4s k - 4p k x 4s, 
- 4p ( while we de all valence configurations 
among the 4s. 4p. and 3d orbital! Including these ad* 
ditional terms tend! to decrease the unponance of 
the duperuon terms and we see a repulsive curve at 
the C AS SC F level but obtain a D t of I 4 kcal mole 
at the l( F level We suspect that Dai MC SCF poten- 
tial curve is too attractive, although his conclusion 
that the experimental results may be reinterpreted in 
temu of a weak bonding model remains valid. 


21 January 191) 

of D f ■ 0.046 eV for 3 1* and 0 061 eV for 1 1* or 
3 X" . This mult is in agreement with the weak bind- 
ing model proposed by Das 

Adi now ledge me nt 

We thank BjOrn Rous for helpful discussions and 
thank Bjftm Roos and Per Siegbahn for use of their 
programs 
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4. Conclusions 


We hive studied the *1. , 1 S* 3 1*, and '-k u 
states of Sc> using C AS SCF CT(SDl wavefunctions in 
a large ga issian basis set We find that the states aris- 
ing out of the *D(4s*3d' )♦ 2 D(4s- 3d* » limit show 
only weak bonding ( *0 06 eV at R * b 0d { ,|. At this 
distance 3d electrons are only weakly coupled and 
the bonding appears to be largely of the van dei Waals 
type arising ojt of the 4s : - 4p* near degeneracy . 

The state which anses out of the ’D(4t*3d* 1 

♦ 4 F 1 4s • 4p* 3d * i limit is found to be st rongly bound 
with respect to -D ♦ 1 1 (4s 1 4p * 3d* M % 0 8 eV at R 
% 7.0 i The bonding in this state is a result of 
interaction between the 4s orbital of the 4s 1 3d* 4p* 
confutation and sp hy brids arising out of the 4s* 3d 1 
and 4s 1 3d' 4p' states States arising out of ; D(4s* 3d 1 1 

♦ 4 F (4s 1 3d* » are predicted to be bound by $ 0.3 eV 
less with respect to *D ♦ 4 F (4s 1 3d* i While this 
atomic asy mptote is too high ( 1 43 eV’l in Sc 2 for 
this state to drop below the 2 D(4s 2 3d * I ♦ 2 D(4s 2 3d * ) 
asymptote, for V, the corresponding asymptote is at 
0.25 eV leading to states from this mixed asymptote 
being the lowest state W'DHG report a significantly 
shorter R t 14.86 a 0 ) and somewhat deeper D t ( 1 12 
eV) for a - Ij state However, our calculations indicate 
similar results for 5 I u and 5 \ . We suspect that the dif- 
ferences between our results and those of W’DHG are at 
least partly due to basis set deficiencies in the calcula- 
tions of W’DHG 

Nt have also studied the binding energy of the 3 1* 
state using the ICF method This leads to the estimate 


References 


1 1 ) C Dai. Chem Phys Utters 86 ( 1982) 482 
|2| C Wood M Moran. I H HUliet and M F. Curst, Fara- 
day S> mp 14(1980) I $9 

1 3 ) J Dro wart and R f Horug.J Phys Chem. 61 U9S7) 
6801 

|4] B uu and A D McLean. J Chem Phys 72(1980) 

3418 

| S ] B Roos and (W Bauschbcher Jr unpublished results 
|6 j A.J H Rachteri.J Gum Phys 58 (1973)4452. 

|7) PJ Hay . J Chem Phy> 66(19*7)4377 
1*1 * Aimlof. MOLECL LE a gausaian integrals program 
|9| RC Ruffenetti, BIGOMOL1, Program No 32*. Quan- 
tum Chemistry Program Exchange, Indiana I'niveruty, 
Bloomington < 1977). 

1 10] P I M Sieghahn. b Roos and C W Baus.hbcher Jr.. 

SWEDEN, a vectorued SCI MC SCI duett Cl 
|II)CW biutchbcher Jr and b H Lengsfteld III Noname, 
a general SCI MC SO open ended Cl program 
|I2)CE Moore, At omit Energy U*els. Cucular 467 
(National bureau of Standards. Washington. 1949). 

1 1 3) C» Biuwhb.her Jr , S P WaJth and H Partridge, J 
Chem Phys 76 (1982) 1003 
1 14 1 S P WakhandCM Bauschbcher Jr . Chem Phys 
Utters 86 < I982i 66. 

CW bauathbeher Jr and S P Walch.J Chem Phys 76 
(1982)4560 

1 1 5 ] B H Botch. TH Dunning h and J E Harruon, 

J Chem Phys 75 (1981) 3466 

1 16 ] SP WaJch unpublished result! 

1 17] P.J Hay and T H Dunning Jr , J Chem Phys 64 <19761 
5077. 

T.H Dunning Jr.SJ* W’alch and A F Wagner. Theoretic- 
al Studies of Selected Reactions in the Hydrogen Oxygen 
System, in Potential energy surfaces and dynamics ealeu- 
lauons.ed DC Truhlar (Plenum Press, New York, 1961), 
and references therein 


ORIGINAL K 

OK POOR O 


1 1 T 


Theoretical evidence for multiple onfr-feiectron 3 d bonding 
in a first row transition metal dimer: The J’ state of Sc 2 

Stephen P Welch* 1 

htlyaiomu j Rtuath Inititult.*' Mountain Fieu. Californio 9404 J 

Charles W Bauechlicher, Jr. 

S4S.4 Amn Ht-teorrh Ctnirr Mo/Jtn Ft eld. California 940 J 5 
(Received 3 Jjai 1913 accepted 26 Jul) 1913) 


Recently, the Sc, molecule had been the subject of 
considerable theoretical study. 1- * The interest in this 
system arises because prev.ous theoretical studies have 
found only weakly bound states for Sc, in contrast to 
mass spectrometric experiments which had been inter* 
preted to indicate strong bonding (D, = 26 1 5 kcal/mol). * 
Previous studies of Sc, b\ Das J and by Walch and 
Bauschlicher 3 indicated that states arising from the 
, D* , D atomic limit were closely spaced and only weakly 
bound with large R, values (~8.0o,). Reinterpreting the 
mass spectrometric experiments u ,ng a larger molec* 
ular degeneracy and larger R, (the original analysis used 
R, = 5. Ma ( , ui # s=230 cm*’, andg=5), Das concluded 
that the mass spectrometric results were consistent with 
weak binding. However, more recently matrix isolatior 
studies have indicated a bound *>' state of Sc,. 1 This 
state appears to be bound with i rg *ct to the ♦ *D 
atomic Limit, but clearly cannot arise from this asymp- 
tote. Walch and Bauschlicher 3 found a ‘a state which 
was bound by - 0. 8 *V with respect to the } V 
+ *F (4s* Ap' 3d') atomic limit, but unbound with respect 
to 3 0 + i D. States arising out of the l D + *F( 4s 1 3 d*) 
atomic limit were no* studied in detail, however, given 
the experimental evidence for a bound *1 state of Sc, 
with a id population of - 3. 0, together with a known 
computational bias toward the 4 F(4s' 4 p 1 3d 1 ) state rela- 
tive to the *F(4s' 3d*) state, ‘ a reinvestigation of the 
*Z* state arising from *D * *F(is l 3d 1 ) was carried out. 


In this communication we report extensive CASSCF/ 
CI(SD) calculations which show that this *1* tcate is 
bound with respect to the *0 * 1 D atomic limit. We find 
that the bonding In this state involves three one-electron 
id bonds and believe this to be strong theoretical evi- 
dence of multiple 3d bonding in a first row transition 
metal dimer. 


The *1* state had previously been suggested as the 
grcund state of Sc. by Wood. Moran, Hillier, and 
Guest 1 however this assignment was not convincing be- 
cause correcting for their asymptotic error the *1^ state 
was not bound with respect to the lowest limit. 


The basis set was the (14sllp6d)' / [8s6p4d] basis set 
described in Ref. 3. An extended basis (14sl2p6d3/) / 
(0*7p4rf2/] was also used. This basis has three 4 /> func- 
tions’ optimized for the V(4s' 4p' 3d 1 ) state of Sc and two 
if functions obtained as a (21) contraction of a three 
term GTO fit to a Slater if with exponent 1 6. The ex- 
tended basis gave an energy improvement of 0.003 39/i 
- 0 092 eV for the CASSCF of the tj sUte at R = 5. Oa 0 . 

As in the previous calculations. 3 the CASSCF space 
consisted of the orbitals derived from the atomic 4s, 
ip, and 3d orbitals. The dominant configuration (r, 

* 0. 85 for CASSCF at R * 5. 0a,) fo, the state is 

4sa}3d<Tj4s«ri3dr'„3diTi. . (1) 

Constraints were placed on the orbital occupancies; 
here the o b'oek had four electrons, the » block had two 
electrons, and the 6 block had zero electrons in all con- 
figurations. The reference configurations for the Cl(SD) 
were selected as those configurations (ten) with Cl coef- 
ficients >0. 05 in the CASSCF wave function near R t . 

This leads to 8TJ reference in the Cl(SD) wave .'unction 
for the *1* state at R = 5. 0a„ ar compared to 8tf? refer- 
ence for the large R asymptote, the 3 E* state at R 
= 50.0a,. The valence populations derived from the 
CASSC F/CI(SD) calculation are 4s = 2. 60. 4p«0.83, 
and 3d - 2. 57. 

Figure 1 shows the calculated CASSCF/CI(SD) poten- 
tial curve for the *1’ state together with the states which 
were previously computed. 3 The *i; state is most 
closely related to the 3 EJ state which has the configura- 
tion 

4*0*4soJldti,3dir‘, . (2) 

As discussed elsewhere, 3 for the states derived from 
the : D * asymptote, the 3d overlap is small and in the 
small S regionone expects one-electron bonds which 
vary with distance like S to be nore favorable than two- 
electron bonds which vary’ with distance like S J . Thus 
the 3 i; state which involves two one-electron t bonds is 
more favorable than the ‘ 1 \ state which has a single 


359C 


J C*WT> Phyi 79(7). 1 Oct 1983 


0021 9606 83 193590^32*02 10 


t 1983 American nrntute ot Pnv*'« 


Letter* to ttie Editor 


3591 




44 . 



FIG. 1. manual entity curve* for the low -lying states of 
Sc ; from CASSCF CHSD) calculations. The locations of the 
: D*'D, : D **F(is' ip ] Xt), and : D * i F{is' Xl 3 i asymptotes are 
Indicated. 


two-electron o bond. Promoting an electron of the 
state from the 4sc, antibonding orbital to the 3 do, orbital 
leads to the 5 I‘ state which has an additioi al one-elec- 
tron 3 dr bond and reduced 4s repulsion due to removal 
of an electron from tbe 4sa, orbital. Given these fac- 
tors. It is not surprist.u’ that the *1^ state has a smaller 


R, and is more stable than the ’l* state 

A Dunham analysis of a parabolic fit to the poten- 
tial leads to R ,* 5. 2^a^ and *, m 1 84 cm* 1 . The binding 
energ> is estimated to be 0. 44 eV (10. 1 kcal 'mol) in- 
cluding (i) a small differential Davidson s correction 
(-0.03 eV). (II) a correction for the error in the asymp- 
totic separation (♦0.26 eV). and (ill) a correction for 
the energv improvement with the extended basis (♦ 0. 09 
eV). The calculated R, and a, values are very close to 
the assumed values. In the bmit of no 3 ft interactions 
*F * ’/) leads to a molecular degeneracy of 280 which 
would lead to an - 16 kcal mol decrease 1 in the experi- 
mental estimate of D, and exrc'lrnt agreement with ex- 
periment. However, the actua' degeneracy will be much 
lower than this, possibly approaching b. and thus the 
calculated l), is smaller than c.irrent experimental esti- 
mates by at least 10 kcal/mol. The calculated w, is 
slightly smaller than the current estimate’ of 238.9 
cm’ 1 which is consistent with the underestimation of D # . 
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The results of extended active space CASSCF calculations are reported 
for Ti^, » an d * Molecular orbitals derived from atomic 4p and 3d' 
are found to be very important leading to improved binding energies for Ti 2 
and V 2 and to a bound curve for Cr 2 • The Cr 2 calculated spectroscopic 
parameters are (experimental values in parenthesis) R e ■ 1.78 A (1.68 A). 

■ 383 cm"^ (480 cm ) and D ■ 0.71 eV (1.56 eV). 
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Recently the transition metal dimers have been of considerable experi- 
mental [1-5] and theoretical [6-10] interest both because of the diverse 
nature of the bonding in the transition metal dimers and because of the 
difficulty in calculating accurate potential curves for those dimers which 
exhibit multiple two-electron 3d bonding (especially Cr 2 )• To date the 
most extensive study of V 2 and has been carried out by Walch et al. [6]. 
These calculations were CASSCF calculations with the 4s and 3d transition 


metal orbitals in the active space. This choice of the active space is 
. . 3 _ 1 4 

appropriate for the states considered, V 2 E anu Cr 2 Z , since they arise 

from the 4s*3d n + 4s* 3d n atomic asymptotes. The calculated CASSCF curve for 


V2 gave good R e and W e values and also these authors correctly predicted 
the ground state of V 2 to be Z prior to experimental determination [4]. 

O 


However, the calculated binding energy for V 2 wa6 only 0.33 eV (relative 
to 4s^3d^+ 4s^3d^) compared to the current experimental value of 1.8 eV 


[4] and C^ did not give a bound curve although a shoulder was observed near 
the region of the experimental R e . 

The usual approach to correcting those defects is the use of config- 
uration interaction (Cl). However, the brute force application of Cl 
techniques to V 2 and Cr 2 is not practical due to the large size of the 
Cl expansions encountered. For example, Walch et al. [6] estimate ^57 
million configurations are needed to describe the potential curve for the Cr 9 
molecule, whereas the current computational capability is only about one 
million configurations. "'Walch et al. were able to carry out a Cl calcula- 
tion for the 1 Z + 6tate of Ti 2 which exhibits a triple two-electron 3d 
bond ( a, n x , 7r v f and used this result to estimate the binding energies of 


V 2 and Cr 2 assuming a triple bonded picture for V 2 , i.e., that the 3d: 


orbita Is are nonbonding. The pre liminary estimate [ 6 ] was an increase in 
binding energy of ^2.0 eV which has now been refined to 1.57 eV leading to 
predicted D e 's of -vl.7 eV for V 2 (relative to 46^3d^ + 4s^3d^) and 0.2eV 
for Cr 2 (relative to 4s*3d"’ + 4s^3d^). This model iB probably realistic 
for V 2 where the 3d 6 orbitals are singly occupied leading to only weak 
one-electron bonds but we now believe it is not reasonable for C^ where 
the 3d 5 orbitals are double occupied (vide infra). 

In order to get around these difficulties we have used a CASSCF 
approach to estimate the dominant correlation effects missing in the CASSCF 
calculations with only 4s and 3d in the active space. For the transition 
metal atoms we find that the most important additional terms are derived 
from 3d' and 4p and these are also found to be the most important addi- 
tional molecular correlation terms. Here 3d' i6 a tight diffuse correla- 
ting orbital for the 3d. Table 1 shows the efrect of including these 


■j. 
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additional correlation terns for the Ti atom. Concentrating on the descrip- 
tion of the 4s 2 3d 2 -» 4s*3d ^ excitation energy we see that SCF underestimates 
the excitation energy. Including the 4p orbital lowers 4s2 relative to 
to 46*and leads to an overestimation of the excitation energy, while inclu- 
sion of 3d' lowers 3d^ relative to 3d? leading to a calculated separation 
that is close to the Cl value when both 4p and 3d' are included. The 
remaining 0.2 eV discrepancy with experiment is due to core-valence 
correlation involving the 3p shell [11]. Thus inclusion of 4p and 3d' 
enables an accurate description of the relative energies of atomic states 
involving different numbers of 4s and 3d electrons. This improved atomic 
description has been shown to be important for describing molecular 
potential curves [12]. An analagous argument for the importance of 3d* has 
been made by Dunning, Botch, and Harrison [13]. 

To the extent that charge transfer terms are important for the transi- 
tion metal dimers, we also expect 4p and 3d' to be important since 4p 
improves the description of terms arising from 

3d n + 4s 2 3d n (1) 

while 3d' improves the description of terms arising from 

4s 1 3d n " 1 + 4s 1 3d n+1 (2) 

Indeed, Goodgame and Goddard [14] have suggested that the difficulty in 
describing the transition metal dimers arises from difficulties in des- 
cribing the ionic terms [especially (2)]. If this explanation is correct 
our calculations include directly the correlation corrections which they 
include by empirical modification of the integrals. 

The dominant configurations [6] for the molecules considered here are: 

Ti 2 o) 

»2 4 6C g3d 0 g3d^ xu 3dTTy U 3d6 X yg3d6 x i_^ (A) 

Cr 2 \ 4s^3d0^d^ u 3dT,2, u 3dJ^ g 3d4_^ ( 5) 

The calculations use the [8s6p4d2f] basis set described previously [6] 

and the MOLECULE-SWEDEN programs [15]. As in the previous calculations the 

occupations are constrained by symmetry. For example for T^ the dominant 

configuration baa four electrons in 6igma, two in tt , and two in tt . 

x y 


For the CASSCF within the As and 3d orbitals we impose the constraints that 
four electrons be distributed among 4sOg, 4sO u , 3d0g, and 3do u , two elec- 
trons be distributed among 3d 7r xu and 3d 7T X g, and two electrons be distributed 
among 3diTy U and 3d flyg. This set of constraints leads to a CASSCF which 
mimics the GVB wavefunction (16). Here we use GVB in the more general sense of 
allowing all possible spin couplings of the As and 3d atomic like orbitals with 
simultaneous orbital optimization. The active space was then augmented by 

addition of 3d* and 4p. Here the additional active orbitals were added 
separately by symmetry blocks. For Ti 2 3do', 3do' plus 4pO, 3d7T^, 3dr^ 
plus ApT^were added in separate calculations. The same occupation con- 
straints were imposed except that in the calculation with 3dTT* plus Ap tt 
the sigma block was constrained to two to four electrons and the 7 ^ block to 
two to four electrons thus allowing the As^ — 4ptf| near degeneracy terms. 
These orbital constraints in addition to keeping the calculations of 
tractable size also lead to a wavefunction which dissociates to Hartree-Fock 
atoms (for the 4s*3d n + 4s*3d n limit considered here). For Cr2 and V 2 
3dC^ , 3dT* and 3d<5^ y were added in separate calculations. For the Ti 2 
calculations the Apa contribution was obtained as the SdC 1 ' plus 4pO contri- 
bution minus the 3dO* contribution and similarly for the 4p7r x contribution. 

The 3d7r', Ap tt, and 3 d6 * contributions were calculated as twice the 3 d 
ApTT x and 3d: xy contributions , respectively. The CASSCF energies and 
energy contributions for each of the added active space orbitals are given 
in Table 2, Table 3, and Table A for Ti 2 , V 2 , and Cr2> respectively. 

Figure 1 shows the calculated potential curves for the Eg state of the 
Ti 2 molecule. The Cl calculation was a multireference singles and doubles 
Cl from nine references leading to ^180,000 CSFs. For comparison to the 
CASSCF curves the Cl energies are shifted to make the asymptotic Cl energy 
equal to the asymptotic SCF energy for the 4s*3d^ + As^3d^ limit (i.e, we 
are comparing binding energies). From Fig. 1 it is clear that both 3d' and 
Ap make important contributions to the binding energy with 3d' somewhat more 
important than Ap. It is also evident from Fig. 1 that the extended CASSCF 
procedure used here obtains a very large percentage of the extra binding 
energy obtained in the Cl. However, estimating Ae 3d' bysumming the 
separately calculated 3d^' and 3 d 7T* components as was done here 
overestimates the effect of 3d' as compared to a calculation in which all 
the 3d' components are included in one calculation. For Ti^ including 3d^', 
3d^ and 3d^ y in one calculation with occupation constraints as above 
increases 3d' by 0.01436h for the molecule at R ■ 3.758Q and by 0.03852h 
for the separated atoms; thus the differential effect of 3d' is reduced from 
' v 32mh to a 8mh. The dominant atomic correlation terms are of the form 

F*ivi :r ; 3d ± - 3d^ * 3dj -3d} (6) 

poor Qu/;Lrr» 
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This result indicates that these terms may not be viewed as "atomic 

correlation" because they are less important for the molecule than for the 
atom# On the other hand, the 4p correlation doe6 not contribute for the 
4s*3d n state of the atom and we do not expect the 4p correlation effect to 
be significantly reduced by a more extensive MCSCF calculation. Given these 
results we must view as fortuitous the good D e for T 12 obtained by including 
3d' and 4p correlation estimates based on calculations which dissociate to 
SCF atoms. However, our intention here is to use this method to extrapolate 
from the Ti 2 case where we are able to do an adequate Cl treatment to the V 2 
and Cr 2 cases where the calculations exceed present computational 
Capabilities. We believe that the correlation effects which are included in 
these calculations will be the most important differential effects between 
molecules in this series. For simplicity in performing the extrapolation we 
use the extended CASCF results directly since the D e obtained for in 
this way is within 0.08eV of the Cl value. 

Table 5 shows the calculated spectroscopic parameters for (obtained 
via a Dunham analysis of the points in Table 2). Here we see that the cal- 
culated is in good agreement with experiment. The calculated D g is smaller 
than the mass spectrometrical ly determined D e but the experimental value 
may be too large as is the case for V 2 where a more recent determination 
via predissociation reduces D e by ^0.7eV. 

Table 6 shows the calculated spectroscopic constants for V 2 . Here the 
effect of 4p was estimated from the Ti 2 calculations. This was done by 
shifting the A E 4p curves for Ti£ toward shorter R by 0.36ao which is 
twice the difference between < r^ s > for Ti 4s*3d^ and V 4s^3d^lll]. From 
Table 6 we see that R e , D e , and w e are all in good agreement with 
experiment for V 2 . In particular the calculated D e is in good agreement 
with the recent experimental value determined from predissociation [4] but 
is considerably smaller than the mass spectrometric value [17]. 

Fig. 2 shows calculated potential curves for Cr2 • Here as for V 2 the 

effect of 4p was estimated from the Ti 2 results but shifted toward shorter 

1 3 

R by 0.64a o which is twice the difference between < r 4 s > for Ti 4s 3d and 
Cr 4s*3d5 [11], The most significant feature of Fig. 2 is the appearance of 
a well in the region where a shoulder had been observed in the CASSCF curve 
[6]. Thu6, Cr2 exhibits rather different behavior from Ti 2 or V 2 in that 
3d' correlation in Ti 2 and V 2 deepened the well but did not significantly 
change R £ or u) e , while for Cr 2 3d' correlation converts a shoulder on the 
curve into a well. This behavior implies considerable R dependence in the 
3d correlation as is evident from Fig. 3 where the 3da , 3d tt , and 3d5 
components of the 3d' correlation are plotted. The 3d 6 component shows an 
especially strong R dependence and it is this term which is mainly respon- 
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»ible for the difference between Cr 2 and V 2 or Ti 2 . In fact from Table 3 
one sees that for V 2 the A E 3d <5 is sma 11 ( ^ 0.0 1 0) and not strongly R 
dependent. Thus a triple bond model based on Ti 2 • a® U6ed in ref 6, is an 
appropriate aero order approximation for V 2 but not for Cr2 • 

From Fig. 2 one sees that there is some enharmonic ity in the calculated 
Cr 2 curve. However, this anharmonicity decreases as the calculation is 
improved and neither the CASSCF 3d’ or the CASSCF + 3d’ * 4p curves show 
any resemblance to the very peculiar potential curve obtained in the recent 
local spin density (LSD) calculations [20]. Also the LSD wavefunction is 
described as antiferromagnetic even near R e whereas our wavefunction shows 
this spin coupling at large R but becomes a multiple 3d bond near R e . 

Table 7 6hows the calculated spectroscopic constraints for Cr 2 • Here 
the best calculation has R e too long by 0.10 A and too small by 97 cm. 
This result is consistent with underestimating the binding energy. However 
it seems probable based on what is known about the binding energy of V 2 
that Cr 2 is bound by considerably less than 1.56eV. The most likely source 
of error in the potential curve is the estimation of the A4p using the 

Ti 2 results. Shifting the R g by 0.64a o based on the relative 4s sizes 
for the Ti and Cr atoms is a somewhat arbitrary assumption although the 
potential is not particularly sensitive to the choice of this parameter 
since very similar spectroscopic parameters are obtained for a shift of 
0.50. It is also possible that the A 4p is slightly larger for than 

for Ti 2 due to the smaller 4s orbital size. However Bagus and Bauschlicher 
report that 4s -*-4p is nearly constant for the atoms varying from 0.76 for 
Sc to 0.85 for Cu [11], Unfortunately, the calculation of A 4p by the 
technique used for Ti 2 would involve a very large CASSCF calculation (which 
we may carry out in the future) since the Cr 2 CASSCF for 3dcr' involves 
^20,000 configurations. 

One further possible source of error in the Cr 2 calculation is core- 
valence correlation [21] involving the 3p shell. This effect would be more 
likely to be important for Cr than for V or Ti since <r 3d > / <r ' : ip > only 
1*39 for Cr while the corresponding ratio is 1.47 for V and 1.56 for Ti 
[11 ) . However, calcu 1 at ions on the 4sl3d^ state of the Cr atom [22] show 
that singles and doubles Cl out of the 4s and 3d electrons leads to a slight 
contraction (6<r> * -0.09a o ) of the 4s and a slight expansion ($<r>*= 0.04aQ ) 
for the 3d, while simultaneous inclusion of single excitations out of the 3p 
has little additional effect. This result makes a significant core-valence 
effect due to the 3p shell seem unlikely. 
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We conclude that inclusion of 3d' and 4p in the valence space leads to 
a very significant improvement in the potential curves for Ti 2 , V 2 , and Cr 2 . 
The calculated R e , D e , and oi e are in excellent agreement with experiment 
for V 2 . For T^ the calculated u e is in excellent agreement with experi- 
ment but the calculated D g is o-O.SeV (with respect to 4s^3d^ + 4s~3d^) 
which i6 smaller than the mass spectrometr ic value of ^1.3eV. However, 
based on previous experience with the mass spectrometric D g for V 2 which 
has been shown to be too large by > 0.68eV, we believe the calculated D fi is 
more reliable. We predict an R fi for T 12 of 1,97 A. To date no experimen- 
tal determination of this quantity exists. For Cr 2 we report the first 

variational calculations with proper dissociation which have the inner 

o 

multiple 3d bonded well. The R e is too long by- 0.10 A and the <n e is too 
small by -»9 7 cm”* but this still represents a significant improvement over 
previous calculations. The calculated D e is 0.71eV which is smaller than 
the mass spectrometric estimate of 1.56eV. Given the large R e and small <u e 
in the Cr 2 potential curve it is quite likely that the true D e is larger 
than the calculated D e although we doubt that it is as large as 1.56eV. 
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Table 1. CASSCF/CI Calculations for the Ti Atom 


4s 2 3d 2 state Asbd 3 state A E (eV) 


SCF 

-.37208 

-.35146 

0.56 

CASSCF 

(4s,4p,3d) 

-.40292 

-.35146 

1.40 

CASSCF 

(4s, 4p, 3d ,3d’) 

-.40756 

-.37072 

1.00 

CASSCF 

(4*#4p,3d,3d f , 4f ) 


-.38366 


Cl a 

-.42234 (-.42649) 

-.38981 (-.39093) 

.89 

exp 



0.81 


a 


Values in parenthesis include Davidson’s correction 


Table 2. Energy Contributions and CASSCF Energies for Tig £ T 



t - 3.50 

3.75 

4.00 

4.25 

CASSCF a 

-.71,240 

-.71684 

-.71051 

-.69821 

AE 3do f 

0.01323 

0.01238 

0.01094 

0.01105 

Ae 4p o 

0.00749 

0.00773 

0.00829 

0.00810 

2* AE 3d tt ' 

/V 

0.01700 

0.01970 

0.02364 

0.02840 

2x AE 4p7r x 

0.01504 

0.01476 

0.01320 

0.01100 


a The CASSCF energies are referenced to -1969.0 hartree 
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Table 3. Energy Contributions and CASSCF Energies for V2 £ 

B 



R - 3.15 

3.25 

3.50 

3.75 

CASSCF a 

-.69183 

-.69666 

-.69407 

-.68236 

A E 3do * 

0.00987 

0.00997 

0.00932 

0. 00769 

2x A E 3d7Tj 

0.01726 

0.01854 

0.02194 

0.02478 

2x AE 3d6 ' 

xy 

0.00870 

0.00924 

0.01058 

O.OU56 


a Tbe CASSCF energies reported here are slightly lover than the values 
reported in ref. 6 which suffered from incomplete orbital optimisation. 
The CASSCF energies are referenced to -1885.0 hartree. 


Table 4. Energy Contributions and CASSCF Energies for Cr 
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Table 5. Cilculitfr# Spe*’troactffU, for Ti 2 



CAS 

* 

CAS + 3d tt 

Exp. 

0 

E g ( A 

1.97 

2.00 

1.99 

1-97 


D e . eV 

0.37 

1.25 

1.86 

XiU 

2 . 9 l a 

w , cm" 1 

436 

423 

442 

438 

407. 9 b 


*Maaa spectrometric value (ref. 17) energy referenced to 4a *3d^ ♦ 4 s* 3<P 
b ref. 19 
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Table 7. Calculated Spectroscopic Constants for Cr 2 


CAS CAS+3d ’ CAS*3d '+4p (est.) Exp. 

O 

R e . A 1.81 1.78 1.68 a 

D e , eV -1.4 0.13 0.71 1.56 b 

W e , cm” 1 308 383 480 c 


aref. 1-3, 

bmass spectrometric value (ref. 18) 
c ref . 19 



Figure Captions 


1 + 

Fig, 1. Calculated potential curves for the state of Ti2. 

1 + 

Fig. 2. Calculated potential curves for the 2 state of Cr ^ 
Fig, 3. Calculated energy contributions for 3d orbitals for Cr 
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Calculated CAS SCF potential cunes are reported for the 3 e" state of V 2 and the state of Cr;. At the CAS SCf 
level the state of >'2 is calculated to be bound ( R e = 1 ,77 A, w e = 593,6 cm" 1 ,D e = 0,33 eV) and to imohe a triple 3d 
bond; while ‘the Crj potential curve is not bound but show s a shoulder near the e.vperimental R e and the wavefunctior, 
shows multiple 3d bonding in this region. 


It is now clear from experimental work that the Cr 2 
molecule has a very short bond length [1], 1 .68 ± 0,01 
A. This short bond has generally been thought to in- 
volve a quintuple 3d bond arising out of the 3do, 3drr, 
and 3d5 orbitals of the Cr 2 4s ] 3d 5 + 4s 1 3d 5 asymptote. 
Given the < r ) for the Cr atomic orbitals (2] <r 4s > = 

1 .94 A, <r 3d > = 0,72 A one sees that at the equilibrium 
bond distance (/? e ) of Cr 2 the 3d orbitals are at a rea- 
sonable distance for bonding but that tliis /? value is 
well inside the optimal bonding radius for the 4s elec- 
trons and they are expected to be non-bonding or anti- 
bonding. 

At larger R values one would expect the spin cou- 
pling to change from predominately 3d~3d bonding 
to 4s-4s bonding, since at large R the 4s orbitals have 
a much larger overlap than the 3d orbitals, The situa- 
tion at large R is expected to be somewhat complex 

4 Supported by NASA grant NCC2-148, 

0 009-2614/83/S 03,00 © Elsevier Science Publishers B.V. 
(North-Holland Physics Publishing Division) 


because the ground state of the Cr atom is high-spin 
6 S arising out of 4s , 3d 5 . As pointed out by Goodgame 
and Goddard [3], forming a 4s— 4s bond with the re- 
maining electrons high-spin-coupled results in the loss 
of a large amount of exchange interaction. As a result 
the lowest molecular state in the large R region is 
found to be an antiferromagnetic singlet state which 
couples the two Cr atoms into a singlet but still pre- 
serves the maximum number of atomic exchange in- 
teractions. 

For Mo 2 Goodgame and Goddard [3] found an 
inner and outer well with the inner well corresponding 
to multiple 3d bonding while the outer well corre- 
sponded to the antiferromagnetic coupling. The per- 
plexing problem is why the Goodgame and Goddard 
Cr 2 calculation did not show a similar 3d bonding solu- 
tion in the inner region, 

In this paper we present complete active space self 
consistent field (CAS SCF) calculations (4J with large 
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gaussian basis sets for the state of V 2 and the 
•Sg state of Ck These calculations show very clear 
evidence for a triple 3d bond in V 2 ; the bonding in* 
volves 3dc and 3drr two-electron bonds with the re- 
maining two electrons involved 3d6 one-electron 
bonds. The bond length here is 1 .77 A compared to 
an experimental value (5] of 1 .757 A and the vibra- 
tional frequency is 593.6 cm” -1 compared to an ex- 
perimental value [6] of 537,5 cm -1 . The calculated 
binding energy at the CAS SCF level is only 0.33 eV 
compared to 2.48 eV for experiment [7] , For Cr^. 
while the wavefunction shows clear evidence of 
quintuple 3d bonding at small R (<1 ,85 A), the cal- 
culated potential curve shows only a marked shoulder 
at %1 .7 A and is siiil repulsive by ss] .4 eV compared 
to current interpretations (8] of the experimental es- 
timate o fD t (1 .0 ± 0.3 eV), At the present time it is 
not computationally feasible to do configuration in- 
teraction (Cl ) calculations for the 3 Sp state of V 2 or 
for the 7 state of Cr 2 at the level of our previous 
calculations on Se 2 |9,10J. An idea of the size of the 
Cl for Cr- may be obtained by noting that a multi- 
reference singles and doubles Cl (MRSD Cl) for Cr 2 
using the basis set with 4f functions and the 3088 
configuration CAS SCF defined below as reference 
configurations leads to 57 million configurations in 
^2h symmetry However, less extensive .MRSD Cl cal- 
culations [II] have been carried out for a 'Sg state 
of Ti 2 which shows triple 3d bonding, 3do and 3d^ 
as in V 2 and Cr 2 , respectively, Here we find an increase 
inD e of * 2,0 eV for Cl as compared to CAS SCF, Ap- 
plying thus correction to the V 2 and Cr 2 CAS SCF re- 
sultsleads to 2.3 and 0.6 eV as our best estimates of the 
D e of V 2 and Cr 2 , respectively. We believe that these 
calculations provide the first convincing theoretical evi- 
dence for multiple two-electron 3d bonding in a first- 
row transition metal dimer. Multiple one-electron 3d 
bonding had been observed previously for Sc 2 by 
Wa),;h anc Bauschlicher [9] . 

The spd basis sets used here were the (14sl I p6d)/ 
[8s6p4d] segmented sets which have been described 
previous!) (10]. These basis sets were augmented 
with two 4f functions, The 4f functions were opti- 
mized based on SDCI calculations for the 4s 1 3d 4 
and 4s ] 3d 5 configurations of the V and Cr atoms 
using a (14si Ip6d3f)/[8s6p4dlf] basis set, The op- 
timum exponents (STO) were 2.7 and 3.2 for V and 
Cr, respectively. The [2 f] contraction is a (21) con- 
traction based on Stewart’s (12] three-term fit of a 
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gaussian f function to a Slater 4f. 

The CAS SCF space for these calculations con- 
sisted of the orbitals derived from the atomic 4s and 
3d orbitals, Constraints were placed on the orbital 
occupancies. For V 2 the o block had four electrons, 
n x and it., had two electrons each, and 6 xy and 
6*: _,,2 had one electron each, For Cr 2 the 6 V) , and 
6*? J,,2 blocks had two electrons each and the other 
occupancies were the same as for V 2 , These con- 
straints basically correspond to a generalized valence 
bond (33] (GVB) treatment of these systems, Here 
we use GVB in the more genera! sense of allowing all 
possible spin couplings of the 4s and 3d atomic like 
orbitals with simultaneous orbital optimization. For 
Cr 2 these orbital constraints reduce the CAS SCF cal- 
culation from 28784 configurations to 3088 configu- 
rations in D 2 h symmetry, The dominant configura- 
tions in the small R region are the SCF configurations; 

4sog3da"3d7rj u 3d7r|. u 3d6 ( |. V g3d6*2 (1) 

for V 2 and 

4so*3d{7g3d7rj u 3dn , ', u 3d6^ l .g3co < j2_ <) ,2g (2) 

for Cr 2 . At R = 3.00 a 0 the percent SCF (C5) in the 
CAS SCF wavefunction is 76 ft for V 2 and 56ft for 
Cr 2 . This clearly indicates the need for an MC SCF 
treatment of these systems. 

Table 1 gives the natural orbital occupancies of the 
CAS SCF wavefunction?. for V, and Cr 2 at R = 3.0 
and 3,5 a Q , In discussing the results in table 1 it is use- 
ful to consider an analogy to the H 2 molecule. Here 
including only left -right correlation leads to a two- 
configuration MS SCF. 

Cj Iso 3 +C 2 Iso 3 . (3) 

Near R e the 1 so 3 (SCF ) configuration dominates but 
at large/?, Cj * C 2 . As discussed elsewhere, in terms 
of the generalized valence bond (GVB) wavefunction 
(13) , the relative weights of the two configurations 
are related to the overlap of the two H(ls) like GVB 
orbitals and with Cj dominant indicative of a 
large overlap <$ f l0 r > whije C j = C 2 implies a very 
small overlap. In interpreting table 1 we car, thus com- 
pare the overlaps of each pair by comparing the natu- 
ral orbital occupancies of the bonding and anti- 
bonding combinations of each atomic like orbital. 
Looking first at Cr 2 , where each possible 3d pair is 
doubly occupied, we see that the above analysis indi- 
cates the 3d overlaps are in the order 


Volume 10 J number 3 


CHI MIC U PUT SICS LITTIRS 
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Tjtlc I 

C AS S( I na'ural orbital ocdipmun lot \ < and ( r ; at A * 3 00 e 0 and J JO« 0 


Stale 

A 

Natural orbital oc.upation *’ 





4so, 

3dOg 

4so u 

3do u 

Jd "*u 


3dt *v* 

?d *x.»u 

v > 

300 

I 93 

1 83 

009 

0.15 

1 90 

0 10 

093 

0.0* 


3.50 

1 93 

1 81 

007 

0 18 

1.85 

0.15 

0.17 

0.13 

Crj 'r; 

300 

1 90 

1 83 

0.11 

0 16 

1 85 

0.15 

1.68 

0.32 


3.50 

1 15 

1 60 

0.15 

0 39 

1.66 

0.34 

1.29 

' l 


*' The 4 m and 3do naiuial orbitals are a mixture of 4i and 3do at A ■ 3 00 « 0 but thou only sbirht amour). > ot mixed character 
at A » 3.5 e 0 


) 


<3drr 3dr> % <3do!3do> ><3d6 : 3d& > . (4) 

Given this analysis one expects the lowest state in V 2 
to arise by distributing two electrons among the 3d6 
orbitals sint.e the 3do and 3drr orbitals are in the 
large overlap region near R ( and the 3d6 orbitals have 
much smaller overlaps in this region Thus it is favor- 
able for the 3do and 3drr electrons to form two-elec- 
tron bonds Furthermore the small overlap of the 
3d6 orbitals makes one-electron bonding more favor- 
able than two-elewtron bonding ( 10) Thus the 
state arising from 3d6|, 3d6j[:_ v : is expected to be 
below the *1* and ’T slates arising out of 3d6*j and 
3d6p_ v i These ?rguments are consistent wi h the 
observation of a -‘I ground state with a short bond 
length for |5| and are also consistent with the 
much weaker bonding in Cr 2 as compared to V, due 
to the necessity of two-electron 3d6 bonding (unfavor- 
able in the small overlap region hn Cr 2 . 

Table : 

Vj ( AS SCI energies 'l' state 


A 

CAS SCI energy 

a.b) 

| 8 v6 p4 d | 

|8v6p4d2f| 

300 

-0.64007 

-0.6 '220 

3 15 


-0.69119 

3.25 

-0.67716 

-0.69621 

3.50 

-068402 

-069390 

3.75 

-0.67798 


4.00 

-0.67051 



*' In au relative to - 1885 0 

b> The a*vmpu tK enerpes (based on atomi. SCF calculi 
ticnv art 4s : 3d' ♦ 4s*3d J - 1885 68586 4«.'3d 4 * 

4s* 3d' - 1885 68076. and 4s*3d 4 ♦ 4»‘3d 4 - l»85.6 1 566 


Tables 2 and 3 give the CAS SCF energies for V, 
and Cr 2 while the results are plotted in figs I and 2. 

A very striking feature of these calculations is the 
large effect of 4f functions in the basis set This ef- 
fect has also been noted by McLean and Liu (14) for 
SCF calculations on Cr : using a large Slater basts set 
Table 4 shows the change in SCF energy at R ■ 3 0a ( , 
as a function of the completeness of the 4i baas set 
The contractions shown in table 4 are based on the 
gaussian fits of Stewart |I2) for an STOof 3.2. The 
largest basis set is the basis set used in the molecular 
calculations. The total eflect of 4t' functions for this 
basis is 1 1 7 eV as compared to I 1 1 eV in the calcu- 
lations of NkLear. and Liu The slightly larger effect 
of 4f functions in our calculations probably reflects 
a slightly poorer spd basis in our calculations A per- 
haps surprising feature o f table 4 is the relatively good 
results with much less complete 4f basis sets In par- 
ticular a single set of 4f primitives obtains 1 .01 eV of 


Table 3 

Crj CAS SCF energies 'i* state 


A 

CAS SCF enerp l bl 

|8s6p4d | 

|8s6p4d2fJ 

3.00 

-0.50554 

-0.53396 

3.25 

-0.53004 

-0.54482 

3.50 

-0.54536 

-0.55196 

3.75 

-0.55917 

-0.56151 

4.00 

-0-47231 



al In au relative to -2086.0 

bl Tht Crj asymptetk energy (based on atomic SCI c*J.ula 
tionsi is -2086 60456 for 4» 1 3d* * 4s 1 »d ! 
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Fi* I Calculated CAS SCF potential curses for the *1* 
state of \ j basic I is the 1 8»6p4d | bauv w hik baa> II is the 
|8s6p4J2f| basis. 


the 1.17 cV effect obtained with the best 4f basis set 
From table 4 one sees that the effect ot 4; functions 
is shghtls smaller at the CAS SCF level than at the 
SCF level 


60- 




R au 


I IK- 2 Calculated C AS SCI potcntul curves fot the 'l* 
state of Ct] ba.M‘ I and basis II are the same size as in ftp I 


Table 4 

Comparison ot 41 bas'i vets fot C tj 't* 


41 batlv 

SC f energ) 

none 

-20M 82 H) 

<D Ml 

• 2085 861*4 

«5» Ml 

-2o8‘ 86190 

mm 

- 2065 86250 

mci 

-208$ 86640 


An interesting feature of fig I is that the effect of 
4f functions is ver> R dependent and shills the /? t 
value of V, to smaller R From fig 2 one also sees 
that the addition of 4t basis functions results tn a dis- 
tinct shoulder on the Cr 2 potential curse in the region 
of the experimental R r While this feature is not evi- 
dent in the calculation without 4f functions, the pear- 
ly linear behavior from 4 00 to 3 .25 a tl is indicative of 
some additional bonding interaction occurring ir. this 
region . 

Comparing to Goodgame and Goddard we note 
that our configuration space contains all the configu- 
rations included in the GVB wave function Our basis 
sel is more extensive than that of Goodgame and 
Goddard which is onls double /eta in the 3d space 
and has no 4f functions As pointed out b> McLean 
and Liu 1 14] a major weakness in the Goodgane and 
Goddatd calculations is the very limited basis set and 
this clearly accounts for missing the shouldet or the 
Cr 2 curse at small R 

Kok and Mall [15] claim to obtain a correct bond 
length for Cr 2 using essentially a perfect pairing GVB 
wasefunction (13) However this result is clearly an 
artifact of the failure of their wavefuncuon to disso- 
ciate properly 

We conclude that the ground states of both V - and 
Cr 2 are characterized by multiple 3d bonding Given 
this, it is not surprising that 4f functions are found to 
be very important in the basis set This new theoreti- 
cal result demonstrates the diverse nature of transition 
metal bonding in that Sc ; has been shown to have 
multiple one-electron 3d bonding |9], V, and Cr- are 
now shown to have multiple two-electron 3d bonding, 
and Cu 2 |I6] has been shown to be predominate^ 

4s -4s bonding with the 3d electrons essentiallv 
atomic like 
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On the nature of the bonding in Cu 2 

Chari** W B«u»chlicher Jf and Stephan P Welch 

Poiyaiomiei <««' xh lutituu * Woo* id in lieu Csli/o'n , « M(MJ 

Per E M Siegbehn 

Imhlltl* 0/ Tktoriitcal Phvua. )« Waidiri J S IIJ 4t SiorkMolm Swede* 
iRrvf 'tO I February I9i: auxpud 26 February 1912 

The ground suit of Cu. * found lo in* from the ‘V4j W '*) ♦ '5|4i W "i Umii and to involve • 4 j Ai 
ugms bond pair Tb« douuMai bond pair correlation* arc left -right and angular with Ihr forme: lengthening 
the bond and the latter Kwitrarur.j the bond, to that at the two-tlactron MCSCF level the A, u slightly longer 
(0 02a,! than at the SCI level Correlation of the W electron! vhortem the bond by 0 14a,. leading to a ftnai 
bond length of 4 lie,, which a 0 lie, longer than eipenment Thu error u of the tame magnitude at twice 
the retail viatic contraction of the 4t orbital of the '5 itate of the Cu atom 10 I Je„i and moat of the remaining 
error in A is thought to be due to thu relativistic contraction 


INTRODUCTION 

Since the ground state of the Cu atom la (4t < 3d") l 
the bonding In Cu, la expected to Involve a 4* -4* bond, 
giving riae to a ground state with all d orbltala fully 
occupied Thta makes Cu, one of the simplest transi- 
tion metal dimers. It la also one of the few tranaltion 
metal dimers for which *iiere ts an accurate determina- 
tion of the bond length (4. 20 bohr), ' For these reasons, 
Cu, serves as an Important benchmark for theoretical 
calculations on transition metal bonding. There have 
been several previous calculations on Cu, (see Rer. 2 
and references therein) and probably the best Is that of 
Peltssier, 1 who used effective core potentials (FCP’s), 
a reasonable valence basis, and a nonvariational Cl 
treatment. One interesting feature of this work was the 
bond shortening al the Cl level when the d electrons are 
correlated. Bagus and BjBrkmar observed this for NIH,* 
and we have observed It for N1H, CoH, and CuH. 4 This 
bond shortening and Its Implications for the nature of the 
transition metal bond are discussed In this article. 

COMPUTATIONAL DETAILS 

The basis set is the lisfybd basis of Wachters, * op- 
timized for the *S state, and augmented with the Hay* 
diffuse d and Wachters’ 1 two diffuse Ap functions, opti- 
mized for i P, scaled by v r i. 25 to make them more suit- 
able for a ground state. This basis Is contracted in a 
segmented manner to ll4.<>ll/>6<f 8«. 6/>4r/]. The 3s com- 
binations of the fb are excluded In some of the calcu- 
lations, an / polarization function is added. This 16 a 
3 GTO fit to an STO with t * 5. 0. The / exponent was 
chosen by extrapolation of the optimum values for the 
Fe and Nt atoms. When the / function is used, only 
seven components are retained. 

All calculations were performed with MGLFCULF — 
noname-Siegbahn Cl, ’ using either the NASA Amts 
7600 or CRAY IS 

Cu ATOM 

In these calculations, symmetry and equivalence re- 
stricted SCF calculations were performed These were 
followed by Cl calculations Including all single and double 
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excitations from all valence configurations with the cor- 
rect symmetry. The 3 d and 4s valence orbitals are In- 
cluded in all calculations, while In the more extensive 
calculations, the 3s and 3 p semicore electrons are also 
included. The calculations In which the 3s and the 3s 
and 3 p are included are denoted as Cl (3/>) and Cl (3s3/<), 
respectively. The results of these calculations are 
summariz'd In Table I. Cl * 3/ denotes that the / functlor 
has been ur«:ontr acted. 

From Table I, one sees that the SCF separation dif- 
fers from : umerlcal HF (NHF) by 0.07 eV, while at 
♦he Cl (0 3 />)■*•/ level, the agreement with experiment 
Is within 0.01 eV. However, when the / function Is un- 
contracted, 3d l0 4s' is stabilized by “0. 1 eV with respect 
to St/Vs 1 , thus worsening the agreement with experi- 
ment. This Is consistent with our previous work, 1 
which Indicated that the 3<r ,1 4>. 1 states of the transition 
metal atoms are stabilized relative to the 3<f4s* states 
upon uncontractlng the / functions Including relativis- 
tic effects stabilizes 4**3 <f by 0.43 eV with respect 
to 4s , 3d ,f '. * Thus, It is not surprising that improvement 
in the calculation leads to 4**3 d* too high with respect 
to 4s , 3</'\ However, we note here that the molecular 
calculations neglect relativistic effects, and a good 
description of the separation is obtained at the 

Cl (3*3 p) + f level. Thus, we chose to use this basis 
set for the molecular calculations. 


TABLE I. Summary of Cu atom calculations. Davidson cor- 
rections are enclosed In parenthesis. 


Calculation 

‘Stf"**') 

’DU/' 4s 1 ) 

A 

SCF 

-1638.613 256 

-1636.796 967 

0.44 

Cl 

- 1639. 0*4 659 

> 1639.015 096 

1.35 (1.30) 

C \*f 

-1639. 166 449 

- 1639. 129 926 

1.59 (1.55) 

0(3/)) */ 

-1639.360 351 

- 2639. 327 964 

1.43 (1.41) 

Cl(3s3 p) *f 

- 1639. 436 39’ 

- 1639. 360 992 

1.51 (1.50) 

C\*V 

-1639. 217 3«- 

->639. 155 872 

1.67 (1.65) 

Cl<3/» 

- 1639.430 515 

- 1639.374116 

1.54 (1.53) 

CUii3p)*V 

-1639.494 622 

- 1639.435 222 

1.62 (1.631 

EXP* 



1.49 

NHF 



0. 37 


*C. E. Moore, Atomic Energy Levels, Natl. Bur. Stand. (VS; 
Clrc. 467, (1949). 
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TABLE n. Surr.msry at bond lei^the lor Cu, calculations (to 
bohrm). Devidsoa corrections art enclosed (a parenthesis. 


Cslculsboe 

Without /’• 

With ft 

•cr 

4.69 

4.53 

Cl 2 electros 

4.42 


2 el MC1CF t«r!-<*L 
2 el MCSC r i«j -<%. vj - 'it 

4. 70 


4. 41 


U frozen) 
«</*e relaxed) 

4. 01 


4 el rVMCfc f me, 4i.4 rj 

4.41 


6 el FVMCaCF (gff,4»4 p) 

4. 57 


6 a Cl (8C F reference) 

4. 54 (4. 55) 

4. 53 (4. 54) 

6 o Cl <MC reference) 

4.54 


U vl Cl 


4.42 (4.40) 

22 el Cl 

4.42 (4.34) 

4. r. (4. J5) 

34 el Cl 


4. 34 (4. 34/ 

Pelletier* SCF no f’» 

4. 54 


NoovertsUoanl Cl 22 el 

4.25 


Lxpertment* 

4.20 



*H.?ferenc» 2. ^Reference ). 


Cu. CALCULATIONS 

The calculated bond lengths lor Cu, are summarized 
tn Table II. The separation between computed points 
wasO. 0 bohr, so the reported H, values are expected 
to be within 0. U2 bohr of the true minima of the model 
used. 

The first calculation is a single configuration SCF 
... 6oJ'3d)7 c*(4S-4S bond) 6o*.(3</)3r 4 (3</)3r 4 (3d)16;(3</) 
Aiejod). 

Using these SCF orbitals, a aeries of Cl calculations 
were carried out including single and double excitations 
from the SCF configuration with progressively more 
electrons being correlated. The Cl's include a 
2 el Cl (7o # ), a 6 c el Cl (6o # , 7o # , and 6o.), a 14 el Cl 
(6o f , lOf, 6o„ 3» f , and 3r,), a 22 el Cl, which includes 
all 3</ and 4s electrons, and a 34 el Cl which Includes the 
3 p as well as the 3 d and 4s electrons. For each point, 
Davidson's correction 19 was computed, added to the 
computed energy, and the minima redetermined. T'.> s«- 
results are reported in parenthesis in Table IL 

In addition to these Cl calculations, a series of 
MCSCr calculations were performed. They Included 
a two-configuration MCSCF, which Included only 7o* 

- 7oJp which allows proper dissociation. A three-con- 
figuration MCSCF which, in addition to the 7oJ-7oi 
configuration, includes 7<rJ - 4v*(4/> ). This MCSCF was 
performed first with the Is, 2s, 3s, 2 p, 3 p, and 3 d 
orbitals frozen at the SCF level, and then with only the 
core fls, 2s, 3s, Ip, and 3p) frozen, L e. , the 3d was 
relaxed at the MCSCF level. The 3 do electrons (6o, 
and 6c,) were also Included with the core frozen at the 
SCF leveL In thl6 calculation, a full valence MCSCF 
ua* performed with the 6 o electrons In the 6o,, 7o /( 

6c. 70. and 4», orbitals. This is 36 CSF's in £,». 

This MCSCF was followed by a 6 el Cl including all 
single and double excitations from all 36 references. 

The final FV'-MCSCF distributes the 6 o electrons Into 
the 6c,, 7c„ 6c„ 7 o„ and all components of the 4 />'*, 

». *• , 8c„ 8c 4tr„ and 4r, 


ORIGINAL PACE 1*: 

DISCUSSION r >F POOR 

At the SCF level, the H, Is more than 0. 3 bohr longer 
than experiment even when / functions are included In 
the bsats. Considering first correlations of the 4* and 
4 p electrons.. Inclusion of the 7 o, orbital lengthens the 
bond. This Is typical of an MCSCF which Includes the 
bonding and antibonding orbitals. The inclusion of 7c* 
-4rJ shortens the bond into agreement with the 2 el Cl. 
only slightly longer than the SCF. The importance ol 
this angular correlation elfect la a result of the near de- 
generacy ol the 4s and 4 p orbitals. Fren though C i Is 
3d"4»' and near degeneracy la not Important lor the Cu 
atom, the simultaneous excitation of 4* - 4/> on each 
atom ia Important for Cu, with a weight approximately 
equal to the Importance of 7c* -7c*. It is Interesting 
to note that this effect was observed (or Lit by Jonsson 
el al. , 11 In Cr, by Coodgame and Goddard, 1 and In our 
lab lor Cs t . ” Goodgame and Goddard described this as 
s van der Waala term, and Jonsson el al. as a near de- 
generacy correlating with the united atom limit. We 
prefer to view this as a near degeneracy effect, allowing 
angular correlation of the o bond. This excitation (7c* 

— 4w*) reduces the Importance of 7c* -7c*. which is a 
bond lengthening excitation, and moves a pair of elec- 
trons into the 4r, uonding orbital. 


Considering now correlation of the 3 d electrons, the 
Inclusion of only the 3 do electrons hat ont> a small ef- 
fect, as is seen by the do el Cl. The inclusion ol the 
3-*s '3d), the 14 el Cl, shows a large bond shortening 
*0. 14 bohr, while Inclusion ol the 16 orbitals *22 el Cl) 
leads to an additional *0.05 bohr shortening, or about 
1/3 of the change of going from 6 tc ’4 el Cl’s. When 
the 3 p electrons are also Included Us el Cl), virtually 
no additional shortening occurs. 

The large bond shortening due to correlation of the 
3r electrons Is thought to result from three effects: 

(I) the correlation resulting from the 4s, -3r interac- 
tion, due to the Importance of the 7o*- 4s* excitation; 

(II) orbital relaxation effects Included ir. the Cl (tide 
inf? a), and (ID reduction of d-d repulsion between the 
evtrtsrs. 

Fffe"‘s (11) and (111) are also present when the 16 or- 
bitals are correlated, while lor 3 do correlation only ef- 
fect (11) la present; thus, the smaller bond contraction 
due to correlation of these orbitals may be rationalized. 

In Cu,, the atoms are 3d'°4s' like, and evaluating the 
8CF energy of the 3d*4s l state using the Sd'Us 1 orbitals 
leads to an excitation energy of 4. 70 or "3.2 eV 
larger than experiment, whereas Individual SCF calcu- 
lations yield a separation of 0. 37 cV. The Inclusion of 
configuration Interaction improves the separation through 
orbital relaxation and correlation, thus lowering the 
3d , 4s t state relative to the 3d ,0 4s' state. Thi6 is oppo- 
site to the effect if separate SCF calculations are car- 
ried out for each state. This reduction of the separa- 
tion allows an increased 3J*4s* interaction or some d-d 
bonding. Wt note h< e that Pelissier’ stated that Cl 
would not mix in more 3d*4s* since Cl increased the 
separation. This wa6 based upon separate orbitals 
for each state ard therefore does not apply In Cu,. 

Similar effects have been observed for NiH. 4 
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It la interesting to note that Pelissier’s* nonvarla- 
tion 22 el Cl lead* to an .1, 0. 14 bohr shorter than our 
result. II the shortening we observe by including 1 func- 
tions at the 22 el C*l level 0 04 bohr) were subtracted 
from his result, he would be in excellent egreement 
with experiment. However, In our lab, we have found 
a useful estimate of the relativistic bond shortening 
is obtained Iron, the relativistic contraction of the va- 
lenceorbilaisobserted by Desclaux. 11 For Cu, this Is 
0. 064 bohr or an expected bond shortening for Cu } of 
*0.13 bohr. If this estimate (sued, our 22 el Cl with 
/'a would be corrected from 4. 35 to 4. 22 bohrs. Con- 
sidering the approximate nature o( this relativistic 
shortening and the use of Davidson's correction to com- 
pute the bond length, this must be somewhat fortuitous. 
However, it does lead us to believe that there are poten- 
tial problems with the FCP of Pellsiier, which dues not 
Include relativistic effects. It Is clear that to resolve 
this problem, 2 FCP's must be developed, one with and 
the other without relatiristic effects included at the same 
level of treatment, l. e. , basis set and correlation. The 
nonrelattvlstic FCP should be able to reproduce our all- 
electron result. 

A final point is the mechanism by which idPis* allows 
d-d interactions. One might assume that a 4s* atom 
would form a repulsive state; however, we have shown 
in ScH u that the 4* -4 p near degeneracy allows a 4.<. 1 
- 4/> l excitation which polarizes the 4sc electrons Into 
the e, orbital, moving them out of the bonding region 
and allowing an increased 3 d interaction. 

CONCLUSION 

The Cu atoms in Cu, are best described as °4s 1 
with : 4«.-45 bond. The Important correlation effects 
for the 4» electrons are the usual bonding to antibond- 
ing excitation and a bonding to 4r ( (4/>) bonding orbital. 

We attribute thl6 to the near degeneracy of the 4s -ip 
orbitals. This effect has been seen In other systems, 
such as Li,, Cr,, and Cs,. We expect this to be a com- 
mon feature of transition metal -transition metal bond6. 


We see a large bond length contraction at the Cl level 
when th.* id electrons are correlated. We attrloute this 
in part to an improved description of the 3k/ , 4s , -3</‘ : 4s' 
separation, as a result of orbital relaxation and cor- 
relation overcoming the strong orbital bias In favor of 
3</ ,e 4«. • . This allows an increased involvement of id*is : 
and the possibility of id-id interactions. Our best 
bondlength is *0.15 bohr longer than experiment and we 
attribute most of this to relativisitc effects, based upon 
atomic calculations of Desclaux. This is J'fferent from 
Pelissier's calculations, which appear to show little 
relativistic effect. The question of relativistic bond 
shortening is left as an unresolved problem. 
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